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VOLUME  II 
1.  INTRODUCTION 

Volume  II  presents  the  detailed  mathematical  analyses 
used  to  obtain  the  results  presented  in  Volume  I*.  The  analytical 
models  to  be  discussed  represent  mathematical  descriptions  of  the 
dome -transducer  acoustic  interaction  problem.  This  problem  involves 
the  simultaneous  solution  of  the  constitutive  equations  for  the 
acoustic  velocity  field  of  a  fluid  and  for  the  displacement  field 
of  a  solid.  A  solution  to  this  problem  is  impossible  to  obtain 
without  assumptions  which  simplify  the  analysis,  while  still 
including  the  dominant  physical  mechanisms.  The  models  to  be 
discussed  represent  an  attempt  to  analyze  the  general  problem  in 
terms  of  a  simpler  framework.  The  framework  is  constructed  on 
assumptions  which  make  the  mathematics  and  computations  tractable 
and  at  the  same  time  yield  information  about  the  physical  system. 

To  a  large  extent,  the  purpose  of  analyzing  the  general  problem 
from  different  sets  of  starting  assumptions  is  to  check  the 
validity  of  the  initial  assumptions  for  a  particular  model.  How¬ 
ever,  this  approach  has  the  advantage  of  isolating  important 
mechanisms  which  are  dominant  in  the  behavior  of  a  system. 

The  difficulties  associated  with  obtaining  a  general 
solution  fall  into  two  categories.  The  first  category  is  insuf¬ 
ficient  knowledge  about  the  properties  of  the  interacting  system 
and  its  boundary  conditions.  The  second  category  is  lack  of  mathe¬ 
matical  techniques  for  treatment  of  the  general  problem.  Certain 
assumptions  are  usually  made  to  obtain  tractable  models  of  radiating 
systems.  These  assumptions  are  as  follows: 

1.  The  system  can  be  described  by  linearized  constitu¬ 
tive  equations  for  both  the  fluid  and  the  solid. 

2.  The  boundaries  of  the  fluid  in  which  the  system  is 
immersed  are  at  infinity. 

*The  nonconcent ^ic  cylindrical  dome  has  been  treated  in  detail  in 
Ref.  [  1.1  and  will  not  be  discussed  here. 
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3.  The  fluid  is  homogeneous  and  inviscid. 

4.  The  transient  effects  are  neglected,  i.e.,  the  time 
dependence  of  all  field  quantities  is  given  by  e~  1 
where  x  is  the  angular  frequency  of  oscillation  and 
t  is  the  time. 

5.  Velocity  control  is  maintained  on  the  transducer 
surface . 

While  the  above  assumptions  yield  a  well-defined  mathe¬ 
matical  boundary  value  problem  for  the  acoustic  field,  it  is  still 
impossible  to  solve  this  problem  for  arbitrary  dome-transducer 
geometry  and  dome  construction.  Rather,  it  is  necessary  to  devise 
a  different  mathematical  model  to  study  each  part  of  the  general 
problem.  The  assumptions  which  are  employed  most  frequently  in 
the  models  are  described  below. 

In  all  the  mathematical  models  except  the  one  devised 
to  study  the  effects  of  dome  structure  (see  Section  7)  the 
following  assumption  is  made. 

6.  The  dome  material  is  homogeneous  and  has  no  supporting 
structure  members. 

It  can  be  shovm  [ 2]  that  on  the  basis  of  Assumptions  1-6 
the  only  geometrical  shapes  for  which  an  eigenfunction  solution 
for  the  dome  motion  can  be  obtained  are  two-dimensional  circular 
cylinders,  three-dimensi.onal  spherical  shapes,  and  infinite  plane 
surfaces.  In  order  to  exter  3  tha  description  to  geometries  other 
than  spheres  or  cylinders,  it  is  necessary  to  make  some  additional 
assumptions  about  the  dome.  Thus  in  addition  to  the  six  stated 
assumptions,  all  models  except  the  thick  cylinder  model  of  Section 
7  are  based  on  an  additional  assumption. 

7.  The  dome  thickness  is  small  in  comparison  to  an 
acoustic  v/avelength. 

Assumption  7  allows  the  use  of  thin  shell  equations  to 
describe  the  response  of  the  dome. 
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To  simplify  the  dome  analysis  even  further,  it  is  pos¬ 
sible  to  use  a  result  obtained  from  the  thick  dome  model  and  from 
previous  analyses  [3].  At  sonar  transmit  frequencies,  the  effect 
of  a  homogeneous  dome  is  due  entirely  to  its  mass,  for  a  limited 
range  of  dome  thicknesses.  This  leads  to  the  final  assumption, 
employed  only  in  the  elliptic  dome  model. 

8.  The  dome  stiffness  effects  are  negligible. 

A  further  removal  of  geometrical  restrictions  to  include 
a  larger  class  of  surfaces  implies  that  eigenfunction  expansion 
techniques  are  no  longer  applicable  for  analysis  of  the  pressure 
field.  It  becomes  necessary  to  use  more  general  analytical  tech¬ 
niques.  The  most  promising  mathematical  technique  for  analysis 
of  general  surfaces  is  based  on  an  integral  formulation  of  the 
interaction  problem.  An  integral  formulation  of  the  interaction 
problem,  based  on  Assumptions  1-8,  will  be  given,  but  the  resulting 
algorithm  has  not  yet  been  developed  as  fully  as  those  previously 
mentioned. 

From  the  previous  discussion,  it  is  evident  that  the 
mathematical  models  available  for  the  interaction  problem  become 
less  realistic  in  description  of  the  dome  as  the  geometrical 
restrictions  are  successively  removed.  In  practice,  computational 
techniques  impose  an  additi jnal  constraint  on  the  configurations 
which  can  be  analyzed.  The  analyses  in  terms  of  eigenfunction 
expansions  are  more  difficult  for  complicated  geometries  because 
the  special  functions  involved  are  more  difficult  to  compute. 

When  the  transducer  and  dome  belong  to  different  separable  geometries, 
as  in  the  elliptic  dome  model,  it  is  necessary  to  calculate  two 
different  sets  of  special  functions  and  establish  a  correspondence 
between  them.  The  integral  equation  formulation  is  size-limited 
and,  to  analyze  realistic  size  bodies,  it  is  necessary  to  give  up 
part  of  the  method's  shape  versatility.  Thus,  even  though  mathe¬ 
matical  techniques  are  available,  in  principle,  to  analyze  a  wide 
variety  of  geometrical  configurations,  the  application  of  these 
techniques  is  still  limited. 
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2.  GENERAL  METHODS  -  INTEGRAL  EQUATION  FORMULATION 

The  following  four  sections  give  a  survey  of  the  tech¬ 
niques  developed  for  a  numerical  treatment  of  a  restricted  class 
of  acoustic  radiation  problems. 

2 . 1  INTEGRAL  FORMULATION  OF  THE  DOME -TRANSDUCER  INTERACTION 

The  purpose  of  this  section  is  to  obtain  a  mathematical 
formulation  of  the  dome-tr ansducer  interaction,  subject  to  the 
eight  assumptions  stated  in  the  Introduction. 

The  implication  of  these  assumptions  is  that  the  acoustic 
pressi.e  field  can  be  determined  by  replacing  the  transducer  and 
dome  by  equivalent  source  distributions  on  surfaces  defined  by 
the  boundaries  of  the  transducer  and  dome.  (Because  of  Assumption 
7  the  dome  can  be  replaced  by  one  surface.)  The  pressure  field 
is  characterized  in  terms  of  the  source  strengths  by  the  require¬ 
ment  that  the  field  satisfy  the  inhomogeneous  wave  equation.  The 
physical  boundary  conditions  on  the  field  at  the  surfaces  then 
determine  the  source-strengths  as  solutions  to  a  pair  of  coupled, 
integral  equations. 

The  most  general  distribution  of  sources  on  a  surface 
for  which  the  field  remains  finite  in  the  neighborhood  of  the 
surface  is  a  superposition  of  simple  (monopole)  and  dipole  sources. 
The  pressure  field  resulting  from  such  a  distribution  on  the  trans¬ 
ducer  and  dome  surfaces  can  be  written: 

p(r)  =  1  jCDt(q)  G(rjq)  +  vfc(q)  G(r|q)  '  dS(q) 

st  q 

+  L  -d(^  TrT  G^lq)  +  *d(q)  G(r|q)r  dS (q)  ,  (2.1-1) 
Sd  3 
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where 


g(H  q)  =  7^- 


j.  e 


ik  S r-q| 


k-ql 


(2.1-2) 


In  (2.1-1),  r  denotes  an  arbitrary  point  in  the  fluid  medium,  and 
p(r)  is  the  pressure  at  that  point.  The  subscripts,  d  and  t, 
denote  respectively  quantities  associated  with  the  dome  and  the 
transducer  surfaces.  The  variable,  q,  denotes  an  arbitrary  point 
on  either  surface.  The  operator  applied  to  G(r|q)  denotes 

q 

the  outward  (toward  the  farfield)  normal  derivative  with  respect 
to  the  source  coordinate.  '<! (q)  is  the  strength  of  the  simple 
source  distribution  and  is  the  strength  of  the  dipole  source 
distribution  at  the  point  q  on  the  relevant  surface.  The  choice 
of  the  source  function  G  in  (2.1-2)  is  based  on  the  condition  that 
the  pressure  field,  p,  satisfy  the  Sommerfeld  radiation  condition 
at  infinity.  (The  time  dependence  is  assumed  to  be  given  by  the 
factor  e~  a  and  this  factor  is  suppressed  in  all  equations.) 

The  wave  number,  k,  in  (2.1-2)  is  the  ratio  of  the  angular  time- 
harmonic  frequency  and  the  sonic  speed  of  the  fluid. 

In  (2.1-1),  the  pressure  field  exterior  to  the  trans¬ 
ducer  is  expressed  in  terms  of  four  unknown  source  distributions 
which  must  be  determined  from  the  given  boundary  conditions.  In 
general  both  the  field,  p(r) ,  and  its  normal  derivative  experience 
jump  discontinuities  across  the  surfaces  St  and  S^.  The  jump  dis¬ 
continuity  of  a  field  variable,  f(r),  across  an  arbitrary  surface, 
S,  is  defined  by  the  relation: 


[[f]]s 


=  Lim  [f(q  +  nc)  -  f(q  -  nQJ  , 
C-0+ 


(2.1-3) 


where  q  is  an  arbitrary  point  on  S,  and  n  is  the  normal  to  S  at 
q  pointing  into  the  farfield.  By  making  use  of  the  singular 
properties  of  G  it  can  be  shown  that 
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:ipj  JSj_  =  ~t(q)  ,  -'n  *  7p]]s^  =  -  -vt(q)  , 


Llp]]s  =  cd(q)  ,  [[n  *  vp]]s  -  -  vd(q)  , 
d  d 


/  o  i  /  \ 

U.  I-4*) 


and  the  jump  discontinuity  across  any  other  closed  surface  is 
identically  zero.  The  relations  (2.1-4)  hold  at  every  point  on 
the  surfaces  for  which  the  principal  curvatures  are  continuous. 

Use  of  Assumption  7  and  the  fact  that  the  normal  particle  velocity, 
v  ,  of  the  fluid  is  equal  to  that  of  the  dome  implies  that 

[[v]:isd  =  [["  ‘  vp]]r-  =  0  ’  (2.1-5) 

where  is  the  angular  frequency  of  the  time-harmonic  oscilla¬ 
tion,  and  o  is  the  fluid  density.  Thus  from  (2.1-4)  and  (2.1-5), 


•d 


=  0  . 


(2.1-6) 


Assumption  8  and  Newton's  law  applied  to  the  dome  implies, 

[[p33s  =  -iu>a  Nd  ,  (2.1-7) 

d 

where  Nd  is  the  normal  velocity  of  the  dome  and  a  is  the  surface 
density  of  the  dome.  Thus  from  (2.1-4) 


cod  =  -  mo  Nd 


(2.1-8) 


=  iu'c  Nd(q)  . 


(2.1-9) 
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Equation  (2.1-8)  gives  the  physical  interpretation  of  the  dipole 
source  distribution  on  the  surface  of  the  dome. 

It  still  remains  to  determine  the  source  distributions 
on  the  surface  of  the  transducer.  There  is  some  degree  of  arbi¬ 
trariness  in  the  determination  of  these  distributions  because  the 
region  interior  to  the  dome  is  not  a  region  of  physical  interest. 
The  nhvsical  boundary  condition  which  must  be  satisfied  on  St  is 
that  determined  by  Assumption  5,  i.e,, 

nt  *  9p^t  +  ntC)  =  im  Nt(^t)  ’  (2.1-10) 

where  NTt  is  the  specified,  normal  component  of  the  velocity  cf 
the  transducer,  and  qt  is  an  arbitrary  point  on  the  transducer 
surface,  the  surface  curvature  being  continuous  at  that  point. 

This  condition,  (2.1-10),  is  not  sufficient  to  determine 
the  source  distributions  and  v  uniquely,  since  according  to 
(2.1-4)  they  are  determined  from  the  jump  discontinuities  of  the 
pressure  and  the  normal  derivative  at  the  transducer  surface. 

Since  the  solution  for  the  pressure  field  interior  to  the  trans¬ 
ducer  is  not  of  physical  interest,  it  can  be  chosen  to  be  any 
arbitrary  solution  to  the  time-independent  wave  equation.  Once 
this  arbitrary  solution  has  been  specified,  the  solution  in  the 
region  of  physical  interest  is  determined  from  (2.1-10). 

One  of  two  choices  is  usually  made  for  the  arbitrary 
interior  solution.  The  Helmholtz  representation,  used  here,  is 
obtained  by  requiring  the  interior  solution  to  vanish  identically. 
Using  this  condition  and  the  relations  (2.1-3)  and  (2.1-10)  in 
(2.1-4)  one  obtains 

<Pt(qt)  =  Lim  p(q  +  n  5)  ,  (2.1-11) 

c  c  ?-0+  c  c 
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and 


-t(qt)  =  -trr  Nt^t)  > 


(2.1-12) 


where  qt  represents  an  arbitrary  point  on  the  transducer  surface. 
Thus,  in  this  representation,  the  strength  of  the  dipole  source 
distribution  has  the  physical  interpretation  of  the  surface  pres¬ 
sure,  and  the  strength  of  the  simple  source  distribution  is  pro¬ 
portional  to  the  transducer  normal  velocity,  i.e.,  the  normal 
component  of  the  fluid  velocity  at  the  surface  of  the  transducer. 

The  other  convention  which  is  in  common  use  and  which 
is  called  the  “method  of  sources"  is  obtained  by  continuing  the 
solution  for  the  external  pressure  field  into  the  interior  of  the 
transducer  region  with  no  discontinuity  in  value,  i.e.. 


°t 


Thus  from  (2.1-4),  the  above  relation  implies  o>t  =  0.  This  amounts 
to  the  statement  that  the  effect  on  the  external  fluid  of  the 
transducer  is  representable  by  a  distribution  of  simple  (monopole) 
sources.  The  strength,  of  these  sources  does  not  have  a 

direct  physical  interpretation.  Mathematically,  this  distribu¬ 
tion,  from  (2.1-4),  can  be  interpreted  as  the  discontinuity  in  the 
normal  fluid  particle  velocity  which  must  exist  across  the  sur¬ 
face  St  in  order  that  the  pressure  field  satisfy  the  wave  equation 
at  all  points  of  space,  except  St,  be  continuous  at  all  points 
of  space,  and  have  a  specified  normal  derivative  on  one  side  of 
the  surface.  For  rapidly  varying  velocity  distributions,  the 
Helmholtz  ^presentation  seems  preferable  for  computation,  hence 
(2.1-11)  and  (2.1-12)  will  be  employed.  Making  use  of  these 
relations  and  also  (2.1-6)  and  (2.1-8)  one  obtains  the  relation 
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G(r!q) 

q 

- (q)  Nd(q) 


-ixc  Nt(q)  G(r|q)  dS 
~  G(rjq)dS  , 

q 


(2.1-13) 


in  which  the  transducer  surface  pressure,  p(q) ,  and  the  dome 
normal  velocity,  N^(q)  a re  unknown.  These  quantities  are  deter¬ 
mined  from  two  coupled,  integral  equations  obtained  from  (2.1-13) 
One  integral  equation  is  obtained  by  allowing  the  field  point  to 
approach  the  surface,  St,  from  the  farfield  side.  Taking  account 
of  the  fact  that  near  the  singular  point  the  dipole  term  behaves 
like  the  solid  angle  subtended  at  the  field  point  by  the  surface, 
one  obtains  in  the  limit 


p(q«)  = 


•Gt) 

~T~~ 


p(q)  37P  G(qtlq)  -i»P  Nt(q)  G(qtjq)'dS 


Nd(q)  G(qt|q)dS 

q 


(2.1-14) 


The  asterisk  on  the  integral  means  that  the  integral  is  the  limit 
of  a  sequence  of  integrals  obtained  by  excluding  the  singular 
point  with  surface  patches  whose  areas  tend  to  zero.  The  remaining 
integral  equation  is  obtained  by  taking  the  limiting  value  of  the 
normal  component  of  the  gradient  of  (2.1-13)  as  the  field  point 
approaches  the  dome  from  either  side. 

Making  use  of  the  relation  (2.1-9)  one  obtains  on  per¬ 
forming  the  limiting  process  that 


N 


i(V  = 


1  < 


92 


10)0 


GWdU)  -  W  Nt(q)  G(qd|q)jdS 

Q  C| 


i  f 

p  - 


Nd(^}  jTn^  G(^)dS  * 


(2.1-15) 


n 

y 
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The  relations  (2.1-14)  and  (2.1-15)  constitute  a  set  of 
two,  coupled,  Fredholm  integral  equations  which  determine  the 
transducer  surface  pressure  and  the  normal  velocity  of  the  dome. 

The  computational  procedures  used  to  solve  these  equations  are 
the  same  as  those  used  to  treat  a  radiation  problem  with  no  dome. 
However,  in  this  case  it  is  necessary  to  subdivide  both  the  trans¬ 
ducer  surface  and  the  dome  surface  to  obtain  a  set  of  algebraic 
equations  which  approximate  the  integral  equations.  From  a  compu¬ 
tational  point  of  view  the  only  details  which  are  new  to  the  inter¬ 
action  problem  are  involved  in  establishing  quadrature  formulae 
for  the  kernels  associated  with  the  dome  normal  velocity.  The 
kernel  multiplying  N^,  in  (2.1-14)  is  not  singular  on  the  dome 
surface  and  hence  it  is  easy  to  treat  numerically.  (The  same  is 
true  for  the  kernels  associated  with  in  (2.1-15).  Moreover 
the  numerical  details  associated  with  the  kernels  defined  on  St 
have  been  discussed  at  length  in  the  literature.)  On  the  other 
hand,  the  kernel  involving  in  (2.1-15)  needs  to  be  defined  more 

precisely  in  order  that  it  have  meaning,  and  this  is  the  remaining 
detail  which  needs  to  be  discussed  in  this  section. 

The  last  term  on  the  right  hand  side  of  Eq.  (2.1-15)  is 
of  the  form 


.2 

1  Cq)  =  |  Kq')  37^'G(5lV)  dS(q')  (2.1-16) 


where  S  is  a  closed  surface  with  a  continuously  turning  tangent 
plane  at  every  point  and  q,  q'  are  points  on  this  surface.  The 
function,  <|>  (q)  is  a  continuous  and  differentiable  function  of  the 
surface  coordinates  and  G  is  defined  in  (2.1-2).  In  the  context 
used  above,  the  integral  on  the  right  side  of  (2.1-16)  is  defined 
to  be: 
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Uq)  Lin  -n  •  /  ^  v(q')  ~f  G(q  +  nC|q')  dS(q')  ;  ,  (2.1-17) 

where  n  is  the  normal  to  S  at  q.  To  analyze  this  limit  it  is 
convenient  to  introduce  the  source  function  for  Laplace's  equation: 


G0(r|q) 


r  -  q| 


(2.1-18) 


Consider  the  function: 


Ffr)  -  ”  §V  dS(q')  , 


(2.1-19) 


where  S  is  the  same  surface  which  defined  I(q).  Mathematically, 
F(r)  is  %tt  times  tne  solid  angle  subtended  at  r  by  the  surface 
S,  and  from  the  divergence  theorem  it  follows  that 


1,  r  inside  b 

F(r)  =  <  ky  r  on  S 

0,  r  outside  S  . 


(2.1-20) 


Thus  from  (2.1-20)  it  follows  that 


vF(q  +  nC)  =0,  C  t  0 


(2.1-21) 


This  means  that  adding  any  multiple  of  F(r)  to  the  right  side  of 
(2.1-17)  will  not  affect  the  limit,  hence  on  multiplying  (2.1-19) 
by  -  *Cq)  and  adding  the  result  to  (2.1-17)  one  obtains 


I(q)  =  Lin,  n-v{s[*(q')|^<;(5+S{|q')  -  ♦  (q)^0(^+5£  |  qO]  dS) 
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or 


I(q)  -  *  (q)  Lim  j  j|^[G(q  +  nClq#)  -  GQ(q  +  nClV)3  dS 

^>  ~*  w  o 


2 

+  C-0  ^s[^(V)  ‘  *(^)3  SnBn,G(^  +  »C|q#)  dS(q')j  . 

(2.1-22) 


To  evaluate  the  first  limit  on  the  right  side  of  (2.1-22)  consider 
the  vector  function 

d(r)  =  7  j _  afr  [G^iq)  -  G0^|q)]  dS  »  (2.1-23) 

S  q 

where  r  is  not  on  S.  The  method  of  evaluation  of  the  right  side 
of  (2.1-23)  depends  on  whether  r  is  exterior  or  interior  to  S 
although  the  final  result  is  independent  of  that  fact.  Assume  r 
is  the  interior  to  S  and  consider  the  region  of  space  enclosed 
by  S  and  the  surfaces  S  of  a  small  sphere  centered  at  that  point 
r.  In  this  region 

V'2  [G(r| r/)  -  Gc(rl r')D  =  -  k2  G(r| r')  .  (2.1-24) 

Integrating  (2.1-24)  over  the  region  bounded  by  S  and  T>,  employing 
the  divergence  theorem,  and  showing  that  the  contribution  over 
S  vanishes  cS  the  sphere  excluding  r  tends  to  zero,  one  obtains 

-  ji-  [  G  (r|  q)  -  G0(r|  q)]  dS  -  -  k2  f  Gfr|r')  dvfr')  . 

S  q  V 

(2.1-25) 
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Thus  from  (2.1-23) 

J(r)  =  -  k2v  v  G  (r|  r')  dV(r')  . 


(2.1-26) 


Again,  surround  the  point  r  by  a  small  sphere  centered  at  a  fixed 
point  rQ  and  break  up  the  region  of  integration  into  two  parts: 
the  volume  V  included  in  the  small  sphere  of  radius  "a"  centered 
at  rQ  and  the  remaining  volume  V-V.  The  integral  over  V  can  be 
performed  analytically  to  yield 

v  _  G(rjr')  dV(r')  =  (ka)2  h^ka)  j,(k|r-r0!) 

for  |  r-rQ|  <  a;  where  and  h-^  are,  respectively,  spherical 
Bessel  functions  of  the  first  and  third  kind. 

For  the  remaining  region,  the  gradient  operation  can 
be  performed  under  the  integral  sign.  Moreover,  because  of  the 
form  of  G,  the  gradient  with  respect  to  field  variables  is  the 
negative  of  the  gradient  with  respect  to  the  source  variables. 
This  fact  and  the  divergence  theorem  allows  the  volume  integral 
to  be  replaced  by  surface  integrals  over  S  and  1>.  Again  the  con¬ 
tribution  over  S  vanishes  as  the  radius  of  the  sphere  tends  to 
zero  and  the  limit  of  the  analytic  result  above  tends  to  zero 
with  ''a"  when  the  limit  is  performed,  keeping  |  r-rQ|  <  a.  Thus 
the  result  is 


(r-rQ) 


r-r 


J(r)  =  k2  6  n(q)  G(r|q)  dS(q)  (2.1-27) 

The  same  result  is  obtained  in  a  more  straightforward  fashion  for 
r  outside  S.  Thus  from  (2.1-23)  and  (2.1-29) 
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n  •  [G(q  4-  nC*q')  -  GQ(q  +  nC|q')3  dS/ 

w 

=  k2  Lim  n(q)  •  n(q')  G(q  +  n§|q')  dS(q')  , 
f— c  '  S 

=  k2  *  J  n  •  n'  G(qlq')  dS(q#)  . 

J  S 


(2.1-28) 


The  asterisk  on  the  last  integral  in  (2.1-28)  has  the  same  inter¬ 
pretation  as  mentioned  previously.  The  equality  between  the  two 
latter  integrals  in  (2.1-28)  is  based  on  the  result  [4]  that  the 
two  limiting  processes  give  identical  results  for  a  surface  dis¬ 
tribution  of  simple  sources  which  is  piecewise  continuous.  The 
latter  integral  on  the  right  side  of  (2.1-28)  shows  that  the  limit 
of  the  integral  on  the  left  exists  for  all  surfaces  S  having  con¬ 
tinuous  curvature  [4]  . 

To  evaluate  the  second  limit  on  the  right  side  of 
(2.1-22)  it  is  convenient  to  introduce  the  expression 


jlL 

d  nS  n' 


— ♦  -  — ♦  j  ) 

G(r|  r  )  =  k  cos  y  cos  y 


G  -  sW  <K>  (2 . 1-29) 


where 


R 


y 


COS  V 


n  •  R 
R 


cos  v' 


The  first  term  on  the  right  hand  side  of  (2.1-29)  is  proportional 
to  a  simple  source  distribution,  hence  from  classical  results  the 
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portion  of  the  limit  corresponding  to  this  term  exists.  Moreover 
the  multiplicative  facto*'  of  the  mixed  partial  in  the  second  term 
is  analytic  in  R.  Thus  an  examination  of  the  second  limit  reduces 
to  the  consideration  of  an  expression  of  the  form 

H(r)  =  n  •  r"  [v(q')  -  £(q)]  Go(r|q')  dS  ,  (2.1-30) 

S  q 

where  Gq  is  defined  in  (2.1-18). 

According  to  classical  results,  if  £(q')  is  continuous 
the  jump  discontinuity  of  H  across  S  is  zero.  However,  it  is 
necessary  that  £  be  twice  differentiable  with  respect  to  the  sur¬ 
face  coordinates  in  order  that  the  limit  exist  on  S.  These  mathe¬ 
matical  restrictions  on  0  will  always  be  satisfied  by  the  source 
density  on  a  dome  which  is  being  driven  by  a  pressure  field  and 
thus  one  can  conclude  that  both  limits  in  (2.1-28)  have  meaning, 
implying  that  the  mathematical  description  of  the  physical  prob¬ 
lem  is  at  least  meaningful.  On  the  other  hand,  the  computation 
of  the  limiting  value  of  H  in  (2.1-30)  is  not  as  straightforward 
because  this  limit  is  not,  in  general,  equal  to  the  limit  of  a 
sequence  of  surface  integrals  which  exclude  the  singular  point. 

For  the  conventional  numerical  approximation  scheme  which  is 
discussed  in  the  literature,  the  knowledge  that  both  limits  exist 
in  (2.1-22)  is  sufficient  to  construct  an  algorithm  for  solution. 
Namely,  if  the  subdivision  of  the  surfaces  is  made  fine  enough 
so  that,  to  the  required  degree  of  approximation,  the  variation 
of  the  source  strength  is  negligible  over  the  individual  surface 
elements,  then  the  contribution  to  the  integral  (2.1-30)  over  the 
surface  element  for  which  the  limit  is  difficult  to  perform  is 
identically  zero  for  all  values  of  r,  hence  its  limit  is  also 
zero. 

Thus  it  has  been  shown  that  the  dome  interaction  prob¬ 
lem,  subject  to  the  above  assumptions,  can  be  described  by  a  set 
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of  couple^  integral  equations  which  can  be  solved  by  the  same 
techniques  used  to  treat  bare  radiators. 

2 . 2  HELMHOLTZ  KERNELS  FOR  QUASI -SEPARABLE  GEOMETRIES 

In  Section  2.1  it  was  shown  that  the  integral  equations 
for  the  source  distributions  involve  kernels  which  are  derived 
from  a  simple  source  function  by  differentiation  with  respect  to 
a  coordinate  normal  to  the  given  surface.  A  numerical  solution 
to  these  integral  equations  for  general  surfaces  involves  a  large 
amount  of  computation  because  of  their  two-dimensional  character. 
The  purpose  of  this  section  is  to  discuss  a  restricted  class  of 
radiating  surfaces  for  which  the  two-dimensional  integral  equation 
can  be  decomposed  into  a  set  of  equivalent  one- dimensional  integral 
equations . 

It  is  known  from  differential  geometry,  that  the  prop¬ 
erties  of  a  surface  are  determined  by  specifying  its  first  and 
second  ground  forms: 

2 

dq-dq  =  ^  ha2  (dua)2  ,  (2.2-1) 

a=l 

2n ' dq  -  ••  l  ha\  (dua)2  .  (2.2-2) 

a  =  l 


The  relations  (2.2-1)  and  (2.2-2)  are  the  ground  forms  written  in 
the  principal  coordinate  system,  and  the  K  denote  the  principal 
curvatures.  The  presence  of  the  negative  sign  in  (2.2-2)  arises 
from  the  convention  that  the  curvature  is  considered  positive  if 
the  surface  bends  away  from  the  normal.  Thus  the  curvature  of  a 
sphere  is  a  positive  constant.  A  general  surface  is  defined  by 
specifying  the  dependence  of  the  scale  factors,  h  ,  and  the  curva¬ 
tures,  K  ,  as  functions  of  the  surface  coordinates.  From  the 
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point  of  view  of  (2.2-1)  and  (2.2-2),  the  most  simple  surfaces  to 
analyze  are  those  for  which  the  scale  factors  and  the  curvatures 
are  independent  of  the  surface  coordinates.  It  turns  out  that 
the  only  surfaces  for  which  these  conditions  are  true  are  for 
infinite  plane  surfaces,  and  for  these  surfaces  the  method  of 
inversion  can  be  used  to  obtain  a  general  solution  of  the  wave 
equation . 

The  next  class  of  surfaces  in  order  of  complexity  are 
those  for  which  the  ground  forms  are  independent  of  one  surface 
coordinate,  say  the  U2  coordinate.  It  will  be  shown  that  the  only 
classes  of  surfaces  for  which  this  is  true  are  either  (a)  infinite 
cylinders  with  arbitrary  cross  section,  or  (b)  axisymmetric  bodies 
of  revolution.  For  surfaces  of  these  types,  the  dependence  of 
the  source  kernels  on  the  symmetry  coordinate  of  the  field  point, 
U2,  and  the  source  point,  u^,  appears  only  in  the  combination 
(u2  -  u p ,  and  it  is  this  fact  which  allows  the  decomposition 
mentioned  previously. 

2.2.1  Quasi-Separable  Surfaces 

In  the  discussion,  the  orthonormal  vectors  tangent  to 
the  principal  coordinate  lines  are  denoted  by  e  ,  the  normal  to 
the  surface  by  n,  and  differentiation  with  respect  to  surface 
coordinates  by  a  subscript  preceded  by  a  comma.  In  terms  of 
this  notation,  the  mathematical  definition  of  the  surfaces  under 
consideration  is  given  by 


=  0 


(2.2-3) 


In  view  of  relations  (2.2-3)  and  (2.2-1)  there  is  no  loss  in 
generality  in  setting  h^  to  unity.  This  gives  u^  the  interpre¬ 
tation  of  being  the  arc  length  along  the  surface  curves  of  constant 
values  of  U2.  Making  use  of  (2.2-3),  the  fact  that  the  u^  are 
principal  coordinates,  and  the  fact  that  ^ ,  n  form  an  ortho¬ 
normal  triad,  one  obtains  the  partial  differential  equations: 
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This  equation  gives  the  interpretation  of  C  and  3;  i.e.,  viewed 
as  space  curves,  the  u?  curves  have  curvature  C/b^  and  a  normal 
which  makes  an  angle  --6  with  e ^  in  the  plane  of  and  n.  The 
binormal  to  these  curves  in  a  right-handed  e-,,  n  system  is 

b2  "  sin  3  +  n  cos  3  ,  (2.2-8) 

and  the  torsion  is  obtained  as  the  magnitude  of  the  derivative  of 
t^.  Using  (2.2-4a),  (2.2-4e)  and  (2.2-6a),  one  obtains 

^2  2  =  €1  2  s*‘n  “  ‘  ni  2  cos  s  =  0  • 

Since  the  curves  have  vanishing  torsion,  they  are  plane  curves. 
Moreover,  from  the  first  of  (2.2-4b),  the  plane  does  not  rotate 
as  u-^  is  varied.  This  means  that  there  are  two,  fixed,  orthogonal 
unit  vectors,  say  i  and  j,  such  that 

2 ( u l ’  u2^  =  1  cos  a  +  I  sin  a  >  (2.2-9) 

where  a  depends  only  on  U£.  Also  from  (2.2-4)  and  (2.2-6)  one 
obtains 

'2,22  +  °2?2  =  0  *  (2.2-10) 

Thus  from  (2.2-9)  and  (2.2-10)  one  obtains 

a  =  C  U2  +  tt/2  , 

where  choice  of  the  arbitrary  integration  constant  is  made  to 
conform  to  uhe  usual  definition  of  the  spherical  coordinate 
system.  Hence 
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-9  =  -i  sin  C  u?  +  j  cos  C  u. 


(2.2-11) 


Using  (2.2-11)  the  expression  for  the  binormal  to  the  curves 


becomes 


b?  =  ixj  =  k 


(2.2-12) 


Also  the  first  of  relations  (2.2-4a)  and  (2.2-4c)  can  be  written 


-n  _  -  2_n  _  t 

dB  “n  *  ~8  “1 


(2.2-13) 


Using  (2.2-8),  (2.2-12)  and  (2.2-13)  one  obtains 


cos  9  ~  +  ®i  sin  e  *  -k  , 


(2.2-14) 


which  when  solved  yields 


-~l  =  "k  sin  S  4-  X  (u2)  cos  3  ,  (2.2-15) 

■—4 

where, 1  is  an  arbitrary  vector  function  of  U2«  From  (2.2-13) 


n  =  k  cos  8  +  t  (t^)  sin  8 


(2.2-16) 


Relations  (2.2-15),  (2.2-16)  and  the  fact  that  e  ^ ,  z  2,  n  are 
right-handed  orthonormal  vectors  implies 


l (u?)  =  i  cos  C  u7  +  j  sin  C  u. 


(2.2-17) 
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Hence  Che  final  solution  for  the  unit  vectors  takes  the  form 


— *  — * 

n  =  i  sin  3  cos  C  U2  +  j  sin  3  sin  C  U2  +  k  cos  3 

=  i  cos  3  cos  C  U2  +  j  cos  3  sin  C  U2  -  k  sin  3 

c  2  =  -i  sin  C  u,,  +  j  cos  C  U2  , 


(2.2-18) 


where  5  and  C  are  defined  by  (2.2-6).  As  mentioned  previously 
C/t^  is  the  curvature  of  the  U2  curves,  k  is  the  (constant)  binor¬ 
mal  to  the  U2  curves.  The  interpretation  of  3  follows  from  the 
first  equation  of  (2.2-18)  as  being  the  angle  between  the  normal 
to  the  surface  and  the  binormal  to  the  U2  curves. 

If  C  is  set  equal  to  unity  in  (2.2-18),  one  generates 
the  class  of  axisymmetric  bodies  of  revolution  by  interpreting  u^ 
to  be  the  rotational  coordinate.  If  C  is  set  equal  to  zero,  one 
obtains  the  class  of  infinite  cylinders  with  the  unit  vector  j 
defining  the  axis  of  the  cylinder.  Moreover  the  relations  (2.2-5) 
show  that  when  the  surface  invariants  are  independent  of  both 
surface  coordinates,  then  both  C  and  3  must  be  set  equal  to  zero, 
hence  the  relations  (2.2-18)  reduce  to  the  equations  of  a  plane. 

An  interpretation  for  general  values  of  C  is  obtained 
by  noting  that  the  arc  length  along  the  U2  curve  is  given  by 

^2^ul^  =  ^2U2 


and,  from  above,  the  radius  of  curvature  of  the  U2  curve  is  given 
by 


1/R2  =  C/h2  . 

Therefore,  the  last  relation  of  (2.2-18)  can  be  written  as 
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(2.2-19) 


Thus  (2.2-19)  shov;s  that  any  non-zero  value  of  C  gives  the  same 
physical  information  as  the  value  C  =  1.  In  summary,  it  has  been 
shown  that  the  only  surfaces  for  which  (2.2-3)  are  satisfied  are 
axisymmetric  bodies,  where  infinite  cylinders  are  a  special  case 
of  axisymmetric  bodies  for  which  the  symmetry  axis  moves  off  to 
infinity . 

2.2.2  Integral  Equation  for  Axisymmetric  Surfaces 

Having  established  the  fact  that  all  the  surfaces  defined 
by  (2.2-3)  are  axisymmetric,  it  is  convenient  to  introduce  the 
rotational  coordinate  ©  defined  by  the  relation 

o  =  C  u2  .  (2.2-20) 


Moreover  it  will  be  convenient  to  consider  the  fundamental  triad 
as  being  the  set  k,  l ,  s2  which  are  the  unit  vectors  of  a  cylin¬ 
drical  coordinate  system.  Equations  (2.2-15),  (2.2-16),  (2.2-17) 
and  (2.2-20)  determine  the  relationships  between  this  triad  and 
the  surface  triad.  In  terms  of  the  fundamental  triad,  the  position 
vector  to  a  point  on  the  surface  can  be  written 

q(u#,  o')  =  V  z'(u)  +  T' (cp')  p'(u')  . 


Also  an  arbitrary  field  point  in  space  can  be  written 


r=kz+£(o)  p  . 


Hence 


R  =  |r  -  q|  =  [(z  -  z')  2  +  o2  -  2pp  '  cos  (cp  -  cp ')  +  p ' 2]  ^ 


(2.2-21) 
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in  (2.2-21)  the  primed  coordinates  refer  to  the  source  point  and 
the  unprimed  coordinates  refer  to  the  field  point.  It  should  be 
kept  in  mind  that  z7  and  o'  are  parametrically  related. 

Making  use  of  (2.2-16)  and  the  expression  for  the  gra 
dient  operator  in  cylindrical  coordinates  one  obtains: 


=  n'  •  =  cos  0'  |r/  +  sin  B 7  X,  ,  (2.2-22a) 

on  o  op 

~  =  n  .  7  =  cos  ?  +  sin  3  .  (2.2-22b) 


Also,  the  differential  surface  area  is 

dS(q')  =  h^(u7)  p 7  du  dp  =  0(u7)  du7  dp7.  (2.2-22c) 

The  special  forms  of  relations  (2.2-21)  -  (2.2-23)  allow  the  decom¬ 
position  of  the  integral  representation  for  the  pressure  field 
given  in  Section  2.1,  as  will  now  be  shown.  For  simplicity,  the 
case  of  one  radiating  surface  will  be  considered,  but  the  analysis 
holds  for  any  number  of  radiating  surfaces  as  long  as  they  are 
axisymmetric . 

From  Section  2.1,  the  most  general  form  for  the  pressure 
field  due  to  a  distribution  of  simple  and  dipole  sources  on  an 
axisymmetric  surface  is 


P(r)  =  |  du7ft(u')r  |^,  $  G (R)  dcp7+  T”  *  G(R)  dp7' 

where  (2.2-24a) 


2tt 


2rr 


ikR 


G<R>  *  TFR  > 


(2 . 2-24b) 
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and  the  remaining  geometric  quantities  are  defined  in  (2.2-21)  - 
(2.2-23).  In  writing  (2.2-24),  the  special  form  of  (2.2-22a)  has 
been  used  in  order  to  interchange  operations  of  integration  and 
differentiation.  The  source  distributions,  §  and  C ,  in  (2.2-24a) 
are  functions  of  both  u7  and  c 7 ,  and  the  condition  that  they  be 
single-values  imposes  the  restriction  that  they  be  periodic  in  o' 
with  period  2 -.  The  most  general  function  of  the  surface  coordi¬ 
nates  which  satisfies  the  periodicity  stipulation  is  representable 
in  a  Fourier  series: 


/  \  into 


$(u',  o')  =  T  $  (u7)  e 
in 


(2.2-25) 


Using  (2.2-25),  a  corresponding  representation  for  C- ,  and  the  fact 
that  R  is  also  a  periodic  function  of  o' ,  one  obtains  a  Fourier 
series  representation  for  the  pressure  field: 


P(r)  =  S  Pm(z,  P)  e 
m 


(2 . 2-2.6a) 


where 


Pm(2,  r) 


du7  2 (u 7 )  *  (u7)  ~,G  +  *  (u7)  G  ;  , 
7  -  m'  7  on  m  mv  7  m  j  * 


Gm(z,  o  |  z7(u7)  p7(u7))  -  ^  JQ  e 


lma  e 


ikR(o.) 


(2 . 2-26b) 


(2.2-26c) 


R(0  =  r  (z  -  z7)2  +  P2  -  2po  7  cos  a  +  p/2]^  .  (2.2-; 
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The  relations  summarized  in  (2.2-26)  characterize  the 
pressure  field  m  terms  of  a  Fourier  series,  each  coefficient  p 
defined  by  an  integral  relation  in  terms  of  the  coefficients  of 
the  source  distributions.  Due  to  the  form  of  the  normal  component 
of  the  gradient  operator  given  by  (2.2-22b),  the  operation  of 
taking  the  normal  derivative  of  (2.2-26a)  does  not  disturb  the  cc 
dependence.  Also,  it  is  to  be  noted  that  the  kernels  are  inde¬ 
pendent  of  the  azimuth  of  the  field  point.  Therefore,  all  of  the 
statements  made  in  Section  2.1  concerning  the  derivation  of  integral 
equations  to  determine  the  source  distributions  can  be  applied 
directly  to  each  of  the  integral  relations  defined  by  (2.2-26c) 
to  determine  integral  equations  for  the  Fourier  coefficients  of 
the  source  distributions.  Once  they  have  been  determined,  the 
pressure  field  is  obtained  from  (2.2-26a)  and  (2.2-26b).  From  the 
point  of  view  of  a  numerical  solution  to  the  one-dimensional 
integral  equations  derived  from  (2.2-26b),  the  remaining  detail 
to  be  considered  is  a  method  of  calculating  the  source  kernels. 

The  physical  interpretation  of  the  source  kernels  is 
obtained  from  the  defining  relation  (2.2-26c).  The  function 
is  the  pressure  field  at  the  point  r,  due  to  a  source  distribution 
of  constant  amplitude,  %nr ,  and  phase,  m,  on  a  ring  of  radius  p' 
centered  on  the  z-axis  a  distance  z'  above  the  origin.  Except  for 
the  degenerate  case  for  which  the  symmetry  axis  moves  off  to 
infinity,  there  is  no  known  closed  form  expression  for  G  ,  and  one 
must  be  content  with  an  infinite  series  or  an  integral  representa¬ 
tion  for  the  G  .  However  one  would  expect  the  singular  part  of 
the  Gm  to  have  the  same  behavior  as  for  the  degenerate  case;  hence 
it  will  be  treated  first. 

2.2.3  Degenerate  Case:  Infinite  Cylinders 

To  obtain  the  limiting  form  of  the  relations  (2.2-25) 
and  (2.2-26)  as  the  radius  p'  of  the  ring  gets  large  it  is  con¬ 
venient  to  introduce  the  arc-length  s'  defined  by 
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s'  = 

°  y 


(2.2-27) 


and  to  refer  the  positions  of  the  field  point  and  the  source 
point  to  a  rectangular  cartesian  coordinate  system  with  the  ori¬ 
gin  on  the  ring  and  the  y-axis  tangent  to  the  ring.  The  relation¬ 
ships  between  the  original  coordinates  (c 7 ,  0')  and  (d ,  0 )  and  the 
new  coordinates  (x7 ,  y7)  and  (x,  y)  are 


x  +  o'  =  p'  cos  9'  ,  x  +  p*  =  o  cos  S  , 


y'  -  P '  sin  9 '  , 


(2.2-28) 


y  =  p  sin 


Using  (2.2-27)  and  passing  to  the  limit  of  large  o7  shows  that  the 
source  densities  are  represented  by  a  Fourier  integral: 


i?,S- 


5(u',  s')  =  o'  fm(u’)  e  p  4(S,)  -  J  io(u\  o)  elos  da  , 

—00 

m  (2 . 2-29a) 

where  a  =  ^-r  ,  becomes  a  continuous  variable  in  the  limit  of  large 

values  of  o7,  and  o7  §  (u7)  -  $  (u7 ,  a)  .  Also  for  large  values  of 
’  me7 

o '  it  may  be  shown  that 


■i  cos  8  Tz'  +  sin  P  dx‘  ' 


(2 . 2-29b) 


Moreover  on  using  (2.2-27)  and  (2.2-28),  and  passing  to  the  limit, 
one  obtains: 

G^(z,  x|z7(u),  x7(u))  =  Lim  p7Gm 

'  p  -» 


1  "  ir  (y7  -y)  eik[(z-z')2  +  (x-x')*  +  (y-y')2]  dy' 

[(z-z7)2+  (x-x7)2  +  (y-y7)2]^ 
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Introducing  th<_  variables 


X  =  [(z  -  z’)2  +  (x  -  x')^"2  , 

(y1  -  y)  =  X  sinh  t  , 
n  rv\  1  f”  „i  X (a  sinh  t  +  k  cosh  t)  , 

tx>  “  7JFT  J  e  ’ 


»  N  2t  h 


and 

one  obtains 


(2 . 2-29c) 


or 


G0(X)  =  i  H0(Jk2  -  oZ'x) J 


(2.2-29d) 


where  Hq  is  the  zero  order  Hankel  function  of  the  first  kind. 
Also,  the  limiting  form  of  (2.2-26b)  becomes 


Pp(z,  x)  =  J  du'  h1(u')  \J0(u')  GaW  +  ♦0(u')  Ga(X^  ' 

(2. 2-29e) 


The  relations  (2.2-29)  are  the  defining  equations  for 
determining  the  pressure  field  on  an  infinite  cylinder.  The 
source  kernels  Gp  defined  in  (2.2- 29d)  are  functions  of  a  contin¬ 
uous  parameter,  a,  and  of  the  magnitude,  X,  of  the  vector  difference 
between  field  and  source  point  projected  onto  the  plane  perpendi¬ 
cular  to  the  axis  of  the  cylinder.  The  closed  form  expression  for 
Gp  is  easy  to  calculate  by  well-known  polynomial  approximations; 
moreover  the  normal  derivative  of  Gp  also  has  simple,  approximate 
polynomial  representations,  and  these  representations  provide  a 
way  of  treating  the  singular  parts  and  the  regular  parts  of  the 
kernels  separately  [  5]  . 
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2.2.4  Series  Representation  of  Source  Kernels 

Unfortunately,  it  has  not  been  possible  to  find  a 
representation  of  (2.2-26c)  which  has  all  of  these  desirable 
computational  features.  There  are  many  representations  for  the 
G  which  can  be  obtained  by  expanding  the  simple  source  function 
in  terms  of  the  eigenfunctions  of  a  separable  coordinate  system, 
but  all  of  these  expansions  have  the  property  that  the  real  part 
is  slowly  convergent  when  the  field  point  is  close  to  the  source 
point.  The  reason  this  happens  is  that  each  of  the  terms  in  the 
series  is  analytic  when  the  field  point  and  source  are  equal,  and 
hence  the  series  must  diverge  to  represent  the  singularity.  From 
the  preceding  analysis,  one  can  ascertain  that  the  dominant  part 
of  the  singularity  of  the  G  should  be  the  same  as  the  zero  order 
Hankel  function,  and,  presumably,  the  singularity  could  be  "sub¬ 
tracted"  out  of  the  series.  However,  the  residual  series  would 
still  have  a  branch  point  singularity,  and  ordinary  numerical 
quadrature  formulae  applied  to  functions  with  branch  points  are 
not  very  accurate.  In  view  of  these  considerations,  it  seems 
best  to  use  an  expansion  in  terms  of  functions  each  of  which  have 
a  singularity  of  the  same  type  as  G  itself  so  that  the  singular 
part  can  be  treated  in  the  same  way  as  the  degenerate  case  of  an 
infinite  cylinder. 

An  expansion  which  has  the  above  feature  and  whose 
terms  can  be  calculated  from  three  term  recursion  relations  is 
developed  in  the  following  manner. 

Introduce  the  variables: 


”  “  RO7)  =  C  (z  -  z')2  +  (p  +  p')2]^  ,  (2.2-30a) 


R(a)  =  n  cos  y(a)  , 


(2 . 2-30b) 
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and  the  identity 


eikn  cos  v  _  ^  (u  +  n)in  +  n (k-n)  (cos  y)  > 


(2.2-31) 


to  write  (2.2-26c)  in  the  form 


.n/ 2  su  u+n 


(kTl) 


(2 . 2-32a) 


where 


,(u)  _  r (u )  (u  +  n)  '*2tt  ima  CnU  (cos  y(a)) 


da  .  (2 . 2-32b) 


In  (2.2-31),  is  the  ordinary  Bessel  Function  and  is  the 
Gegenbauer  polynomial  of  degree  n.  The  integral  in  (2.2-32b) 
can  be  evaluated  by  elementary  methods  by  using  the  representa¬ 
tion  of  the  Gegenbauer  polynomials  in  terms  of  hype rgeome trie 
functions : 


C^v)  (cos  y)  =  — cos  Y  2fi 


From  (2.2-26d)  and  (2.2-30)  one  obtains 


'i-B,  *+S±i;  U+%;  Sln2  Y) 


(2.2-33) 


2  2 
sin  y  =  v  cos  a/2  , 


(2 . 2-34a) 


where 


V-  Y’ 

■n 


(2 . 2-34b) 
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Substituting  (2.2-33),  (2.2-34a)  into  (2.2-32b),  integrating  term 
by  term,  and  making  use  of  the  identity 


3F  j  eln°  (cos  a/2)2t  da 
o 


one  obtains  the  result 


0  ,  i  <  m 

1  (  11  \  ,  ^ 

3T  4  +  J  1  1  «  . 


/  \  7 (u)  (j+n)  (2u )  (j+X  -1-1)  (-X  ) 

1±>(V)  =  —7 — 11  J.5.\  —  2al  <$» 


mn 


2nl  (u  +  %)m  m! 


^F2(o+Xnfm+l,  nrt-%;  u-hn+%,  2m+l;  v)  , 


(2.2-35) 


where 


,  _  n  -  1 

Xn - T~  > 


(2.2-36) 


and  the  function  on  the  right  side  of  (2.2-35)  is  the  generalized 
hypergeometric  function. 

The  relation  (2.2-35)  can  also  be  written  in  the  equi¬ 
valent  form 


T^')  (v)  = 
n,m  J 


(u)  (u+n)  (2u)  d 

-TiTTniT -  v  W  3F2(u+Xn+1’  *Xn’  ^  v+*’  v>  • 

(2.2-37) 


From  (2.2-36),  X  is  an  integer  for  odd  values  of  n;  hence  the 
hypergeometric  function  in  (2.2-37)  reduces  to  a  finite  polynomial 
Using  this  fact  one  obtains  the  result  that 


r- 
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T 


(") 

2k+l  ,m 


(v) 


=  0 


for  k  <  m  . 


Also  from  (2.2-26c)  and  (2.2-32a)  it  is  evident  that  the  (v) 
for  odd  values  of  n  represent  the  imaginary  part  of  G  ,  which 
is  non-singular.  On  the  other  hand,  the  real  part  of  Gm  is 
singular  at  R  =  0.  From  (2.2-26d),  (2.2-30a)  and  (2.2-34b)  the 
singular  point  corresponds  to  the  value  v  =  1,  and  hence  the  singu¬ 
lar  nature  of  G  must  be  reflected  in  the  T^'^(v)  for  even  values 
m  nm 

of  n  and  v  =  1. 

The  singular  properties  of  the  generalized  hypergeo¬ 
metric  functions  have  not  been  studied  as  completely  as  those 
of  the  ordinary  hypergeometric  functions,  and  thus  it  is  desirable 

to  examine  under  what  conditions  the  above  relations  for  the  ^ 

mn 

reduce  to  the  ordinary  hypergeometric  functions. 

The  fundamental  identity  given  by  (2.2-31)  holds  for 
all  values  of  u  which  are  not  negative  integers,  and  thus  the 
value  of  v  in  the  relations  (2.2-31)  through  (2.2-38)  can  be 
chosen  arbitrarily  subject  to  the  above  restriction.  On  the  other 
hand,  the  generalized  hypergeornetric  function  is  reducible  to 
a  simpler  form  in  two  ways:  (a)  by  choosing  y  such  that  a  numera¬ 
tor  parameter  is  equal  to  a  denominator  parameter,  or  (b)  by 
moving  the  singular  point  to  infinity  to  obtain  a  confluenc  form. 

The  first  approach  seems  more  straightforward,  and  from 
(2.2-37)  it  follows  that  the  proper  choice  of  u  which  is  indepen¬ 
dent  of  n  is  the  limit  as  u  approaches  zero  through  positive 
values.  Making  use  of  the  relation 


Lim 


'■'(”)  (v+n)  (2j) 


n 


i-o 


+ 


T 

n . 


n  =  0 
n  ^  0 


n 


J 


one  obtains 
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Tran(v)  "  2Fl(Xn+1»  "Xn ;  ®fl;  v)  '  (2.2-38) 


where  the  superscript  v  =  o  has  been  deleted,  i.e.,  T  (v)  = 

CO)  nm 

T^' (v)  .  The  hypergeometric  function  on  the  right  hand  side  of 

(2.2-38)  is  reducible  to  a  finite  combination  of  associated 

Legendre  functions,  i.e., 


...  /  \  n  m  /  d  v  m  f  .  1  —  V\  m  /  2  «  in  / .  «  \  ~ 

Tn,m(v>  =  2 < I+>.  )  (->.  )  v  W  -  >  P).  C1'2'')-  i 

’  N  n'm  n'm  n 

(2.2-39) 


however,  this  relation  is  not  useful  for  calculation  except  for 

m  =  0  because  the  apparent  square  root  branch  point  singularity 

does  not  exist  in  the  T  .  For  m  =  0,  (2.2-39)  reduces  to 

ran 


(v)  =  P  (l-2v)  , 
An 


(2.2-40) 


where  is  the  ordinary  Legendre  function.  The  relation  (2.2-40) 
shows  that  T  n  are  calculable  from  three-term  recursion  relations 
which,  in  view  of  (2.2-36)  and  the  known  recursion  relations  on 
Legendre  functions,  "?5o  not  couple  even  and  odd  values  of  n.  The 
functions  T^  q  for  odd  n  are  ordinary  Legendre  polynomials  and 
for  even  n  they  are  Legendre  functions  of  half-integer  order. 

It  is  known  that  the  three-term  recursion  relations  for 
the  Legendre  functions  are  computationally  stable  with  both 
increasing  and  decreasing  values  of  the  order  if  the  magnitude  of 
the  argument  is  less  than  unity,  so  that  the  representation, 


:  e  „  J  (kn) 
n  _•  n  nv 

Go  =  T  1  ~ 

n=o 


'  n-1 

~T 


(l-2v) 


(2.2-41) 
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obtained  by  substituting  (2.2-40)  into  (2.2-32a)  is  suitable  for 
computation.  (From  (2.2~34b),  v  =  0  corresponds  to  the  case  for 
which  either  the  source  point  or  the  field  point  is  at  one  of 
the  poles  of  the  a:<isymmetric  body.  If  one  uses  the  fact  that 
the  Legendre  functions  of  unit  argument  are  unity,  the  relation 
(2.2-41)  can  be  summed  by  specializing  (2.2-31)  for  v  -  o  and 
%  =  o  to  get  the  familiar  result.  As  mentioned,  v  =  1  corresponds 
to  the  singular  point  so  that  the  physical  range  of  the  argument 
of  the  Legendre  functions  is  between  plus  and  minus  unity.)  More¬ 
over,  it  will  be  shown  that  the  singular  part  of  (2.2-41)  can  be 
isolated  completely.  The  convergence  of  (2.2-41)  is  determined 
by  the  argument  of  the  Bessel  functions,  which  are  known  to 
decrease  quite  rapidly  in  magnitude  after  the  order  becomes  , 
roughly,  twice  the  argument.  Moreover,  the  Bessel  functions 
can  be  calculated  by  three-term  recursion  relations  which  are 
stable  in  the  direction  of  decreasing  n. 

To  establish  recursion  relations  satisfied  by  the 

T  for  arbitrary  integral  values  of  m,  '.t  is  convenient  to 
n,m  J  °  7 

rewrite  (2.2-38)  in  the  equivalent  form 


T 

nm 


(v) 


m 


*T  Un-1 


(v) 


(2 . 2-42a) 


where 


(X+l)  (-X)  v™ 

U™(v)  =  - ~R(2m+iy - F(-\+m,  X+m+1;  2nri-l;  v)  ,  (2.2-42b) 

and  the  subscripts  on  the  nypergeometric  function  have  been 
deleted  since  the  discussion  will  be  restricted  to  the  ordinary 
functions  with  two  numerator  parameters  and  one  denominator  param¬ 
eter. 
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By  making  use  of  the  recursive  properties  of  the  hyper¬ 
geometric  functions  the  following  relation  can  be  obtained: 

A1  F(a-1,  b+1 ;  c;  v)  +  A_2  F(a,  b;  c;  v)  +  A3  F(a+1,  b-1;  c;  v)  =  0  , 

(2.2-43) 

where 


A1  =  b(c-a)  (b-a-1)  , 

A2  =  (b-a)f  (2a-c)  (b-a-1)  -  2a(c-a-l)  +  (b-a-fl)  (b-a-l)v]  , 

A3  =  a(c-b) (b-a+1)  . 

Making  use  of  (2.2-42b)  and  (2.2-43)  one  obtains  the  recursion 
relation 

X  0  +1)  2-m2]U  ®+1  -  (2X+l)[tn2+X  (X+l)  (l-2v)]U™  +  (X+l)  (X  2-m2)uJ'_1  =  0, 

(2.2-44) 

The  relation  (2.2-44)  reduces  to  the  recursion  relation  satisfied 
by  the  Legendre  functions  of  argument  (l-2v)  when  m  is  set  equal 
to  zero.  Moreover,  it  has  the  same  behavior  as  the  relation  for 
Legendre  functions  for  large  values  of  X,  and,  therefore,  it  has 
the  same  stability  properties.  Thus  the  relations  (2.2-44), 
(2.2-42),  (2.2-36)  and  (2.2-32a)  (specialized  to  u  =0)  define 
the  source  kernels,  G  ,  by  a  series  whose  terms  are  calculable 
from  three-term  recursion  relations. 

It  still  remains  to  isolate  the  singularity  at  v  =  1. 

To  accomplish  this,  it  is  necessary  to  apply  a  linear  transforma¬ 
tion  on  the  argument  of  the  hypergeonetric  function  in  (2.2-42b), 
The  result  is 
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u™(v)  =  In(l-v)  V‘“(v)  -  W“‘(v)  , 


,m- 


(2.2-45a) 


where 


V®(v)  =  vm  -----  F(-X+m,  X+m+1;  1;  1-v) 


(2 . 2-45b) 


„ .  .  (-X+m)  _  (X+m+l) 

w?1(v)=v1111  — — 1 —  - -n - L  2$  (n+l)-Cr  (-X+mfn) -$  (X+mfn+l)]  (l-v)n 

X  -  (n!)2 


R=0 


(2.2-45c) 


and  v (z)  is  the  Digamma  Function. 

Both  V™(v)  and  W^(v)  are  analytic  in  the  neighborhood  of 
v  =  1  according  to  the  definitions  (2.2-45b  and  c) .  Thus  (2.2-45a) 
shows  that  the  singularity  of  the  U™  is  logarithmic.  For  integral 
values  of  X,  the  U™  are  non-singular,  as  mentioned  previously. 

By  applying  (2.2-43)  to  (2.2-45b),  it  is  easy  to  show 
that  the  V™(v)  satisfy  the  same  recursion  relations  as  the  U™(v) , 
and  since,  according  to  (2.2-45a),  the  W^*(v)  are  a  linear  combina¬ 
tion  of  the  U™(v)  and  the  V™(v)  ,  it  follows  that  they  also  satisfy 
the  same  recursion  relation.  This  result  allows  the  singular 
part  of  the  Gm  to  be  separated  quite  easily  from  the  non-singular 
part,  and  this  feature  will  be  used  later  in  establishing  numeri¬ 
cal  quadrature  formulae. 

Summarizing  the  relevant  formulae,  one  finds  that  the 
Gm  can  be  written  (for  the  case  u  =  0  in  (2.2-32))  as 


Gm(zp!zV>  =••  k*i 

n=o 


n 


(2.2-46) 


where  n  and  v  are  defined  in  (2.2-30a)  and  (2.2-34b)  respectively. 
The  functions,  u!r(v)  ,  are  determined  by  the  recursion  relation 

A 

(2.2-44)  and  initializing  values  calculated  from  (2.2-42b).  In 


representation  of  U™(v) 
(2.2-45a)  where  the  V™ 
relation  (2.2-44)  as 
the  U“‘(v)  ,  and  the  initializing  values  are  obtained  from  (2.2-45b) 
and  (2 . 2-45c)  . 

Finally,  it  is  necessary  to  consider  the  effect  of 
applying  the  normal  derivative  operators  defined  in  (2.2-22)  to 
the  representation  of  given  by  (2.2-46).  In  terms  of  the 
independent  variables  -  and  v: 
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the  neighborhood  of  the  singularity,  the 
for  non-integral  values  of  X  is  given  by 
and  the  W^(v)  satisfy  the  same  recursion 


3n' 


3  b_V 

^  t)  on' 


(2 . 2-47a) 


3  _  2_2.  JL  +  iz  JL 

3n  3n  ? Tj  dn  Sv 


(2. 2-4 7b) 


The  operator  ~  applied  to  (2.2-46)  gives  a  linear  combination  of 
Bessel  functions,  i.e., 


-jn(k-o. 


so  that  the  sum  corresponding  to  applying  this  operator  to  (2.2-46) 
involves  only  Bessel  functions  and  U-functions.  However,  in  order 
to  evaluate  the  normal  derivatives  it  is  necessary  to  find  an 
expression  for  the  derivative  of  the  U-functions  as  a  linear 
combination  of  functions  of  the  same  type.  An  examination  of 
the  normal  derivatives  of  v  shows  that  they  can  be  written  in 
the  form: 


n  V 

■’  n ' 


sm  g 

— n — 


'  -  2^'-z)cos  !'  +..V.-Pjl8in  e"jv(1_v)  _ 

(.z'-z)  2  +  (o'-o)2  J 


(d ’-p) 


(2.2-49a) 
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5 v  _  fsin_£ 

3n  L  c 


oCz-z^cos  g  +  (p-p')sin  g 

(z-z')2  +  <P-P’)2 


v(l-v)  . 


(2.2-4%) 


The  quantities  in  the  brackets  are  finite  everywhere  the  curvature 
is  continuous  on  the  surface  of  the  radiator,  so  that  for  a  smooth 
body  it  is  only  necessary  to  examine  the  effect  of  the  operator 
v(l-v)^  on  the  U-functions.  (The  motivation  for  specializing 
consideration  to  this  operator  rather  than  the  derivative  operator 
itself  is  seen  from  an  examination  of  the  special  form  of  the 
kernel  Gq  given  by  (2.2-41).  The  derivative  of  a  Legendre  func¬ 
tion  has  a  higher  order  singularity  than  a  Legendre  function 
itself  and  hence  cannot  appear  in  the  physical  description.  On 
the  other  hand,  the  operator  v(l-v)—  applied  to  a  Legendre  func¬ 
tion  gives  a  linear  combination  of  Legendre  functions.  Finally, 
it  may  be  noted  that  the  right  side  of  (2.2-49)  vanishes  for  plane 
or  spherical  radiators  when  the  field  and  source  points  are  on 
the  surface.) 

Making  use  of  (2.2-42b)  and  the  properties  of  the  hyper¬ 
geometric  functions  it  may  be  shown  that: 

dUm  f  1 

V(1_V)^  =  ^-J[m2+X2(l-2v)]Uj1  -  (X2-m2)U“  A  .  (2.2-50) 


The  relations  (2.2-47)  through  (2.2-50)  when  used  in  conjunction 
with  (2.2-46)  define  the  series  representation  for  the  dipole 
source  kernels  which  appear  in  relations  (2.2-26). 

2.2.5  Suirana  ry 

In  summary,  it  has  been  shown  that  the  class  of  axi- 
symmetric  radiators  are  characterized  by  surface  ground  forms 
which  are  independent  of  one  of  the  principal  coordinates.  The 
class  of  infinite  cylinders  is  a  degenerate  case  of  axisymmetric 
bodies  for  which  the  symmetry  axis  is  at  infinity.  For  axisymmetric 
radiators,  the  two-dimensional  integral  representation  of  the  pres¬ 
sure  field  can  be  decomposed  into  a  set  of  one-dimensional 
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integral  representations,  and  the  source  densities  associated 
with  a  particular  axial  mode  can  be  obtained  by  solving  a  one¬ 
dimensional  integral  equation.  The  resultant  three-dimensional 
pressure  field  is  then  obtained  by  a  superposition  of  the  contri¬ 
butions  from  all  the  modes.  The  source  kernels  associated  with  a 
particular  mode  have  a  logarithmic  singularity  and  a  prescription 
for  calculating  these  kernels  has  been  given. 

2 . 3  NUMERICAL  SOLUTION  OF  ONE-DIMENSIONAL  INTEGRAL  EQUATIONS 

In  Section  2.2,  it  was  shown  that  the  Helmholtz  integral 
representation  for  the  pressure  field  due  to  an  axisymmetric 
radiator  can  be  decomposed  into  a  set  of  modal  representations, 
defined  by  (2.2-26),  which  relate  the  modal  amplitudes  of  the 
pressure  field  to  the  corresponding  amplitudes  of  the  source 
densities.  The  amplitude  of  a  particular  mode  of  the  source 
density  is  determined  from  a  one-dimensional  integral  equation 
which  is  derived  from  (2.2-26b)  by  allowing  the  field  point  to 
approach  the  surface  in  the  manner  discussed  in  Section  2.1.  The 
general  form  of  these  integral  equations  can  be  written  as 

b 

cp(u)  =  t(u')  K(u,u')  du'  -  f(u)  ,  (2.3-1) 

a 


where  u  is  the  variable  characterizing  points  on  the  generating 
curve  of  the  radiator.  The  kernel,  K,  and  the  inhomogeneous 
function,  f(u),  are  determined  by  the  methods  outlined  in  Sections 
2.1  and  2.2. 

The  problem  to  be  considered  here  is  that  of  replacing 
the  integral  equation  by  a  set  of  algebraic  equations  which  are 
obtained  by  a  generalization  of  the  usual  methods.  The  basic 
idea  involved  in  this  numerical  approach  is  that  the  infinite 
set  of  relations  summarized  by  (2.3-1)  is  replaced  by  a  finite 
number  of  relations,  i.e.,  the  relation  (2.3-1)  is  required  to 
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hold  only  at  discrete  points  (particular  values  of  u)  on  the  inter¬ 
val  of  integration.  The  fundamental  question  to  be  answered  is: 
given  a  numerical  solution  to  (2.3-1)  which  is  obtained  by 
requiring  (2.3-1)  to  hold  at  N  discrete  points  (referred  to  here 
as  mesh  points) ;  how  large  must  N  be  in  order  that  the  approxi¬ 
mate  solution  is  a  good  representation  of  the  exact  solution? 

It  will  be  shown  in  the  sequel  that  the  approximate  solution 
satisfies  an  integral  equation  of  the  same  type  as  (2.3-1)  except 
that  the  inhomogeneous  term  differs  from  that  in  (2.3-1)  by  a 
function  which  vanishes  at  all  the  mesh  points.  The  details  of 
an  approximation  scheme  based  on  piecewise  continuous  interpolating 
functions  will  be  presented  and  qualitative  criteria  for  the 
choice  of  N  will  be  given. 

Let  the  mesh  points  be  designated  by  a  discrete  index 
ranging  from  1  to  N  and  labeled  in  such  a  way  that 

u  . ,  >u  ,  m  =  1  ...  N  -  1  .  (2.3-2) 

m+l  m 

Moreover,  let  ;p(u)  be  an  approximate  solution  to  (2.3-1)  in  the 
sense  that  it  satisfies  (2.3-1)  at  discrete  mesh  points,  i.e., 

rb  - 

°m  =  !  CD<U')  K(VU,)  du’  ■  m  »  (2.3-3) 

3 

where 

cd  s  m(u  ) ,  f  =  f  (u  ) . 
m  v  m  ’  m  v  nr 

Define 


6(u)  =  co(u)  -  cd (u)  ,  (2.3-4) 

where  cp  is  the  exact  solution  to  (2.3-1).  Then  on  making  use  of 
(2.3-1)  and  (2.3-3)  one  finds  that 
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,b 

c (u  )  -  * (u')  K ( u  , u 1 )  =  0  . 


m 


m 


(2.3-5) 


Finally,  if  one  defines  Che  inhomogeneous  function: 


•~(u)  =  1  (u  )  -  5(u')  K  (u  ,  u  1 )  du'  , 

"  a 


(2.3-6) 


one  finds  on  using  (2.3-4)  and  (2.3-i)  that 

.  b 

r ( u )  =  ~(u')  K(u,u')  du'  -  [f(u)  -  e(u)J  .  (2.3-7) 

a 

Moreover,  from  (2.3-5)  and  (2.3-6), 

-(u  )  =  0  .  (2.3-8) 

nr 

The  relations  (2.3-7)  and  (2.3-8)  show  the  result  stated  previously, 
namely,  that  the  approximate  solution  satisfies  the  integral  equa¬ 
tion  with  a  modified  inhomogeneous  term.  Moreover,  (2.3-5)  and 
(2.3-6)  show  that  the  difference  between  the  approximate  and  exact 
solutions  satisfies  an  integral  equation  of  the  same  type,  having 
an  inhomogeneous  term  which  vanishes  at  the  mesh  points. 

The  nature  of  the  difference  function,  5  (u)  ,  depends  on 
the  particular  type  of  approximation  technique  which  is  used  to 
obtain  the  solution  r  (u)  .  The  discussion  in  this  section  will 
be  limited  to  approximation  schemes  which  are  based  on  the  use  of 
piecewise  continuous  functions.  Specifically,  the  set  of  piece- 
wise  functions  which  are  to  be  discussed  consist  of  polynomials 
which  are  constructed  to  interpolate  the  functional  values  at 
the  mesh  points. 

The  mathematical  definition  of  the  interpolating  poly¬ 
nomials  to  be  used  is  the  following:  Let  u^  be  an  arbitrary 
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interior  mesh  point  and  define  the  midpoints  of  the  segments 
adjacent  to  this  mesh  point  by 


aj  =  (uj-l  +  uj)/2  >  (2 . 3-9a) 

bj  -  (uj+1  +  Uj)/2  .  (2. 3-9b) 


Denote  by  P.(u)  the  interpolating  polynomial  of  degree  2K  which 
is  defined  over  the  interval  <  u  <  b^  .  Then, 


2?  .  0.  (u) 

P  (U)  =  2.  't  (u-uj  Q!  (u  ) 

t=j  -K  1  2  1 


where 

j+K 

Qj  (u)  =  n  (u  -  uk) 

k=j-K 


(2.3-10a) 


(2. 3 -10b) 


and  a  prime  on  denotes  differentiation  with  respect  to  u. 

The  relation  (2.3-10)  is  the  Lagrangian  interpolation  formula  of 
numerical  analysis,  and  (2.3-10)  is  a  restatement  of  the  fact 
that 

Pj(uk)  =  £k  ,  k  =  j-K  ...  j+K  .  (2.3-11) 

It  is  to  be  noted  that  even  though  the  interpolating  function 
reproduces  the  functional  values  at  2K+1  mesh  points,  its  interval 
of  definition  extends  only  a  half  mesh  point  spacing  on  either 
side  of  the  central  mesh  point  u^ .  Also,  the  definition  (2.3-10) 
holds  only  for  those  mesh  points  for  which  K  <  j  <  N-K.  The 
remaining  mesh  points  must  be  treated  somewhat  differently. 
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Alternatively ,  one  can  define  (2.3-3)  to  hold  also  for  the 
remaining  mesh  points  by  extending  the  mesh  past  the  physical 
range  of  the  variable  u  and  defining  the  values  of  cr  in  the  non¬ 
physical  range  in  some  convenient  manner.  This  approach  amounts 
to  making  some  statements  about  the  derivatives  of  the  function 
v  at  the  end  points  of  the  physical  range  of  u.  In  most  physical 
problems  the  end  points  of  the  interval  correspond  to  special 
points  for  which  additional  information  about  the  function  is 
available,  and  hence  can  be  used.  A  particular  example  of  immediate 
interest  is  the  integral  equation  for  the  source  density  on  an 
infinite  cylinder.  For  this  case  u  parameterizes  a  closed  curve, 
and  the  source  density  must  be  single  valued.  This  implies  that 
the  source  density  should  be  extended  into  the  non-physical  region 
as  a  periodic  function. 


Generally,  it  will  be  assumed  that  the  non-physical 
values  of  ©  have  been  chosen  in  a  manner  consistent  with  the 
physical  nature  of  the  function  ©  and  the  parameter  u,  and  (2.2-3) 
will  be  taken  to  hold  for  1  s  j  s  N.  Then,  in  terms  of  the  set 


of  piecewise  continuous  functions, 
can  be  written 


the  approximating  function 


N 

"(u)  =  ^  U[(b  -  u)(u  -  a.)]  P  (u)  ,  (2 . 3-12a) 

J=i 


v’here  l’(t)  is  Che  Heaviside  step  function,  i.e., 

'0  t  <  0 

U(t)  =  <;  (2 . 3-12b) 

[_1  t  >  0 


42 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  '8721 


Making  use  of  the  fact  that 

i^u.)  =  ~  ,  (2.3-13) 


and  of  the  properties  of  the  step  function,  it  is  easily  verified 
that 


1. .  .N 


(2.3-14) 


hence  oc(u)  reproduces  the  functional  values  at  all  the  mesh 
points.  It  is  to  be  noted  that  (2.3-14)  holds  for  any  integrable 
function  P^ (u)  which  is  constrained  by  the  single  condition  given 
by  (2.3-13).  A  choice  of  the  (u)  which  is  different  from  that 
of  (2.3-10)  can  only  affect  the  difference  function,  6 (u) ,  defined 
by  (2.3-14). 

Substitution  of  (2.3-12)  into  (2.3-3)  yields 
N  b. 

~m  *  I.  •  J  Vu'>  K<VU’>  du'  -  “m  •  (2.3-15) 

3-1  a3 

The  relation  (2.3-15)  is  valid  for  arbitrary  integrable  functions, 
and  at  this  point  one  could  proceed  to  derive  a  general  set  of 
algebraic  equations  for  the  cp  by  choosing  an  arbitrary  set  of 
functions  Pj (u)  which  are  normalized  according  to  (2.3-13).  The 
choice  of  interpolating  polynomials  is  based  on  the  fact  that  they 
are  relatively  complete  and  their  choice  does  not  unduly  compli¬ 
cate  the  evaluation  of  the  definite  integrals.  Thus  on  using 
(2.3-10)  one  obtains: 
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m 


N  j+K 


j  =  l  t=j-K 


K 

I  mjt 


m  ’ 


(2 . 3-16a) 


where  the  weighting  factors  are  defined  by: 


b. 


K 

m  jt 


r  J  CL(u') 

]  K(u  ,u’)  -7~  r  ■  \  t -r  -  t1  du 1 

Ja.  m  Qj(ut) 


(2 . 3-16b) 


To  reduce  (2.3-16a)  to  a  determinate  set  of  relations,  it  is 
necessary  to  give  explicit  forms  to  the  whose  indices  are  out¬ 
side  the  range  from  1  to  N.  For  physical  problems  of  interest 
here,  the  most  general  defining  relation  takes  the  form 


•-t+l 


Jl, 

\ 

k=l 


<-)  r 


Mk  ®k 


1  =  1...K 


(2 .3-17a) 


N+-C 


N 

V" 

k=l 


A (+)  -  / 

Ik  'k  1 


=  1...K 


(2 . 3-17b) 


Making  use  of  (2.3-17)  in  (2.3-16a)  obtains: 


N 


^m  Z  Fmk  ^k  ~  Fm  ’ 

k=l 


(2 . 3-18a) 


where 


J,(k) 


JK^  K-t+1  N 

Fmk  Z  Kmnk  +  Z  LAik  Z  Kmn ,  -1+1  +  Al Ik  Z  Km,n,-H-N  ’ 

n=J1(k)  1=1  n=l  n=l+N-K 

(2 . 3 -18b) 
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and 


J-^(k)  =  max  (1,  k-K)  ,  (2.3-18c) 

^(k)  =  m*-n  C^+K,  N)  .  (2.3-18d) 

The  solution  of  (2. 3- 18a)  determines  the  approximate  function  at 
the  chosen  mesh  points.  Once  these  values  have  been  obtained, 
the  relations  (2.3-17),  (2.3-12),  and  (2.3-10)  determine  the  solu¬ 
tion  at  every  point  on  the  interval. 

The  next  step  involves  establishing  some  criterion 
to  determine  how  accurately  the  approximate  solution  represents 
the  exact  solution,  cp(u).  Basically,  there  are  three  types  of 
errors  introduced  in  obtaining  a  solution  to  (2.3-18a).  The  first 
type  of  error  results  from  attempting  to  model  a  continuous  func¬ 
tion  by  the  quasi-continuous  representation  given  by  (2.3-12a). 

The  second  type  of  error  is  introduced  by  the  assumption  that  the 
exact  function  cp  is  representable  by  a  polynomial  of  degree  2K 
in  the  interval  aj  <  u  <  b^ .  The  third  type  results  from  numerical 
calculation  of  the  matrix  elements,  F  and  the  solution  of  the 
algebraic  system  (2.3-18a).  An  error  of  the  third  type  can,  in 
principle,  be  minimized  by  efficient  numerical  algorithms,  as 
will  be  discussed  in  the  next  section.  It  is  to  be  noted  that  the 
difference  function,  6 (u) ,  defined  by  (2.3-4)  accounts  only  for 
errors  of  the  first  and  second  type. 

It  may  be  shown  that  errors  of  the  first  type  vanish 
identically  at  the  mesh  points  if  the  kernel  and  the  inhomogeneous 
function  satisfy  certain  restrictions.  To  see  this,  it  is  neces¬ 
sary  to  substitute  (2.3-12a)  into  the  integral  and  compare  the 
result  to  the  function  cp  (u)  ;  i.e.,  define 
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oL(u)  =  ^  " (u ' )  K(u , u 1 )  d  '  -  f(u)  .  (2.3-19) 


Making  use  of  (2.3-15)  and  (2.3-10)  it  follows  that 


P.(u)  =  f  cp(u' )K.  (u,u' )  du'  -  f.(u)  ,  (2 . 3-20a) 

J  ^  o  J  J 


«  Q,(u) 

Kj(u’u'}  =  X  K(ul*u,)  (u-up-QT^ 

l=j-K  1  J  1 


(2 . 3 -20b) 


Vu)=  I 


3  -K 


Q,  (u) 

f(V  (u-Ulj  q:(u-t! 


(2.3-20c) 


Substitute  (2.3-20a)  into  (2.3-12a)  and  subtract  the  result  from 
(2.3-19)  to  obtain: 


[rL(u)-?(u)j  =  j  v(u’)[K(u,u')  -  K(u,u ' ;  ]du'  -  [f(u) 


-  f(u)]  , 

(2 . 3-21a) 


where 


K(u,u')  =  ^  U[ (b.-u)(u-aj)3  Kj(u,u')  , 

j=l 


(2 . 3-21b) 
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The  relations  (2.3-20)  and  (2.3-21)  show  that  if  the  kernel  and 
the  inhomogeneous  function  are  bounded,  then  that  part  of  the 
difference  function,  6 (u) ,  which  is  due  to  an  approximation  of 
the  first  type,  depends  on  the  error  introduced  by  representing 
the  u-dependence  of  the  kernel  and  the  inhomogeneous  function  by 
an  interpolating  polynomial  of  degree  2K  in  the  intervals  a^  <  u  < 
bj .  Moreover,  these  relations  show  that 


^(uj)  -  <5(Uj)  =  0  ,  j  =  1. .  .N  ;  (2.3-22) 


i.e.,  an  error  of  the  first  type  vanishes  at  the  mesh  points  if 
the  kernel  and  inhomogeneous  functions  are  bounded  at  the  mesh 
points.  It  turns  out  that  (2.3-22)  is  also  satisfied  if  the 
kernel  K  is  merely  integrable.  It  is  important  to  note,  however, 
that  the  relation  (2.3-22)  is  not  sufficient  to  insure  that  the 
solution  cp  (u)  is  a  good  approximation  to  the  exact  function  cp  (u) 
even  at  the  mesh  points.  The  reason  this  is  true  can  be  seen  from 
(2.3-21a).  If  the  left  hand  side  of  (2.3-21a)  is  not  a  small 
quantity  at  every  point  on  the  interval  a  £  us  b,  then  the  origi¬ 
nal  assumption  about  the  variation  of  co  has  been  violated  and 
hence  the  solution  for  the  cpm  is  not  valid.  Thus  from  (2.3-20) 
and  (2.3-2?)  it  follows  that  a  solution  obtained  by  the  above 
scheme  can  be  meaningful,  only  if  the  mesh  points  are  spaced  close! 
enough  that  the  interpolating  polynomial  used  is  of  sufficiently 
high  degree  to  represent  the  variations  of  both  the  kernel  and  the 
inhomogeneous  function  over  the  entire  interval.  This  condition 
insures  that  an  error  of  the  first  type  has  a  small  effect  on  the 
function  at  the  mesh  points. 
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It  still  remains  to  determine  the  effect  of  an  error 
of  the  second  type.  In  order  to  do  this,  it  is  necessary  to 
examine  the  remainder  term  associated  with  a  polynomial  approxima¬ 
tion  of  the  type  given  by  (2.3-10).  A  discussion  of  this  problem 
is  given  in  most  texts  on  numerical  analysis  [6]  .  The  result  for  a 
function  which  is  differentiable  can  be  written: 


where 


a>  (u)  =  P.(u)  +  6 . (u)  ,  a.  <  u  <  b.  ,  (2.3-23a) 

J  J  J  J 


(2. 3-23b) 


and  where  P . (u)  is  defined  by  (2.3-10)  with  the 
the  co-p  f  is  a  point  in  the  interval  Uj_^  <  §  < 


cp  ^  replaced  by 
Uj+K’  anC*  t^ie  SuPPr“ 


script  on  <o  denotes  the  order  of  the  derivative. 


The  relation  (2.3-23b)  shows  that  if  the  function  has 
derivatives  up  to  order  2K+1,  then  the  error  in  a  polynomial  inter¬ 
polation  vanishes  at  the  mesh  points.  This  means  that  in  the 
absence  of  any  weaker  constraints,  an  approximation  scheme  using 
interpolating  polynomials  is  limited  to  those  integral  equations 
which  admit  solutions  that  are  differentiable  up  to  order  2K+1. 


An  analysis  of  the  perturbation  effects  of  adding  a 
term  of  the  form  of  (2.3-23b)  to  (2.3-12a)  requires  further  study. 
One  qualitative  statement  which  can  be  made  here  is  based  on  the 
faut  that  the  perturbing  term  given  by  (2.3-23b)  is  a  polynomial 
of  degree  2K+1.  Thus,  one  could  partially  account  for  this  term 
by  increasing  the  degree  of  the  interpolating  function  in  (2.3-10) 
to  2K+2  by  including  K+l  adjacent  mesh  points  on  each  side  of  the 
central  mesh  point.  Then  one  could  repeat  the  same  arguments  used 
to  derive  (2.3-20)  and  (2.3-21)  to  obtain  analogous  expressions, 
except  that,  in  this  case,  the  kernel  and  the  inhomogeneous  function 
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would  be  interpolated  by  polynomials  of  degree  2K+2.  One  would 
expect  that  if  the  kernel  and  the  inhomogeneous  function  are 
reasonably  well  represented  by  an  interpolating  polynomial  of 
degree  2K  then  the  difference  between  the  approximate  solutions 
obtained  using  interpolating  polynomials  of  degree  2K  and  2K+2 
should  be  small.  Hence,  it  would  appear  that  a  mesh  point  distri¬ 
bution  which  is  fine  enough  to  interpolate  the  kernel  and  the 
inhomogeneous  function  with  a  polynomial  of  degree  2K  is  both 
necessary  and  sufficient  to  insure  that  an  approximate  solution 
based  on  an  interpolating  polynomial  of  the  same  degree  is  reason¬ 
ably  accurate. 

When  applied  to  the  integral  equation  for  pressure 
sources  on  a  radiator,  the  above  considerations  imply  that,  for 
velocity  distributions  which  are  rapidly  varying  (e.g.,  discon¬ 
tinuous)  ,  the  "Helmholtz"  formulation  is  preferable  to  the  "source" 
formulation  since  the  inhomogeneous  function  of  the  Helmholtz 
formulation  is  more  slowly  varying  than  that  for  the  source  formu¬ 
lation. 

2.4  EVALUATION  OF  THE  MATRIX  ELEMENTS 

T.n  Section  2.3,  it  was  shown  that  the  algebraic  system 
of  equations  which  define  an  approximate  solution  to  a  one¬ 
dimensional  integral  equation  depend  on  a  set  of  weighting  factors 
of  the  kernel,  defined  by  (2.3-16b).  In  Section  2.2  it  was  shown 
that  the  kernels  associated  with  axisymmetric  radiators  have  a 
logarithmic  singularity,  consequently  ordinary  quadrature  formulae 
are  inaccurate  when  used  to  evaluate  the  weighting  factors  near 
the  singularity. 

It  will  be  shown  in  this  section  that  the  weighting 
factors  are  linear  combinations  of  a  set  of  quantities  called 
moments  of  the  kernel,  and  numerical  quadrature  formulae  which 
account  for  the  singular  nature  of  the  kernel  will  be  given. 
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From  Section  2.3,  the  definition  of  the  weighting 
factors  is  given  by 


l  mn 


m 


m 


dii' 


K(u  u')  ),.“v 

4  ^  "V 


(2.4-1) 


From  (2.3-10),  the  coefficient  of  the  kernel  under  the  integral 
sign  is  a  polynomial  of  degree  2K  in  the  integration  variable, 
u' .  Moreover,  the  relation 


(u* 


u  ) 
m 


m+K 

V 

zli 

n=m-K 


(U  n  '  um>k  V1'5 

^"r  V  W 


(2.4-2) 


is  an  identity  for  05  k  5  2K  and  k  integral.  Thus  on  multiplying 
(2.4-2)  by  K(u^,  u'),  integrating,  and  using  (2.4-1)  one  obtains 

m+K 

Mlmk  ’  Z  <Un  '  um)k  Ktan  -  k  =  0....2K  -  (2-«) 

n=m-K 

where  , 

b 

„  m  , 

M/mk  ’  J  <u'  '  um)  K(“i.u')  du'  •  (2.4-4) 

am 

The  relation  (2.4-4)  is  the  defining  relation  for  the  moment  of 
degree  k  of  the  kernel  centered  around  the  mesh  point  u  .  The 
relation  (2.4-3)  is  a  system  of  (2K+1)  equations  whose  solution 
determines  the  weighting  factors  in  terms  of  the  moments.  A 
closed  form  solution  of  the  system  of  equations  (2.4-3)  is  obtained 
by  expanding  the  right  side  of  (2.4-1)  and  using  (2.4-4): 
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2K 

^mn  -  Zj  ^mnk  ^-tmk  ’ 

k=0 


where 


C 


=  1 

mnk  kTQTry 


dk  V*> 
<X'V 


(2.4-5) 


(2.4-6) 


The  relations  (2.4-4),  (2.4-5)  and  (2.4-6)  show  that  the 
weighting  factors  are  obtained  by  a  product  of  a  matrix  Cmn  which 
depends  only  on  the  mesh  point  distribution,  and  the  moment  matrix 
which  depends  only  on  the  local  character  of  the  kernel  in  the 
neighborhood  of  the  mesh  point  u  . 

An  evaluation  of  the  moments  defined  by  (2.4-4)  should 
account  for  the  singular  nature  of  the  kernel.  In  Section  2.2 
it  has  been  shown  that  the  general  form  of  an  axisymmetric  kernel 
can  be  written: 

K(u„,  u')  =  S^(u')  In  w^(u')  +  R^(u')  ,  (2.4-7) 

where  from  (2.2-26d),  (2.2-30a),  (2.2-34b)  and  (2.2-45a) 


w^(u') 


(z^z')2  +  (pt-o')2 

(zt-z')  2  +  (P^+P  ')2 


(2.4-8) 


The  singular  point  of  the  kernel  occurs  at  u'  =  u^ .  It  is  con¬ 
venient  to  consider  the  functions  and  as  being  defined  by 
two  different  expressions  depending  on  the  value  of  u' .  Thus,  if 
u'  is  not  close  to  u^  ,  is  defined  to  identically  zero  and  R^ 
is  defined  by  applying  (2.2-47a)  to  (2.2-46)  and  using  (2.2-48) 
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through  (2.2-50).  If  u'  is  close  to  u, ,  the  relation  (2.2-45a) 
is  used  to  isolate  the  logarithmic  part,  and  the  and  R,  are 
defined  uniquely  by  this  decomposition.  Both  representations  of 
the  Sj  and  the  R,  defined  by  this  scheme  are  analytic. 

The  degenerate  case  of  infinite  cylinders  can  be 
treated  in  the  same  manner.  For  this  case,  it  can  be  seen  from 
(2.2-29)  and  from  the  form  of  the  Hankel  function  for  small  argu¬ 
ments  that 


w2(u')  =  (k2-o2)  [  (z^-z')  2  +  (xrx')2]/2  .  (2.4-9) 


Using  the  decomposition  defined  by  (2.4-7),  the  moments  can  be 
written 


^•f.mk  ^Imk  +  R4mk  * 


(2.4-10a) 


where 


m  , 

’-Cmk  =  (u'-um)  lnwi(u')  du'  ,  (2.4-10b) 

m 


and 


b 

ro 


(u''um)k  Rl(u')  du'  ' 


(2.4-10c) 


The  quadrature  formulae  which  have  been  used  for  the  evaluation 
of  (2.4-10b)  are  based  on  an  approximate  evaluation  of  the  inte¬ 
grals 


Tk^x)  =  J  u ^  F(u)  ln  w(u) 

o 


(2.4-11) 
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where  F(u)  is  an  arbitrary  cubic  and  w(u)  is  a 
function  of  u  which  vanishes  only  at  the  lower 
gration  interval.  The  approximation  to  can 


slowly  varying 
limit  of  the  inte- 
be  written: 


Tk(x)  =  xk+1[Tkl  F(0)  +  Tk2  F(gkx)  +  Tfc3  F«]  .  (2.4-12) 


The  values  of  6k  and  the  T^n,  n  =  1...3  are  determined  from  the 
fact  that  (2.4-12;  should  hold  for  arbitrary  cubics,  F(u).  It  is 
straightforward  to  show  that 


n  _  Wk+3  "  Wk+2 
k  ‘  Wk+2  -  Vl 


,  _  Wk+2  '  \  Wk+1 

k3  1  -  8k 

7  _  Wk+1  "  Wk+2 

Lk2 


Tkl  _  Wk  "  Tk2  "  Tk3  » 


( . 24-13a) 

(2.4-13b) 

(2.4-13c) 

(2.4-13d) 


where 

Wk  =  — J  In  w(u)  du  •  (2.4-14) 

v  X  o 

The  evaluation  of  the  moments  will  be  discussed  presently. 
First  it  is  convenient  to  point  out  that  under  the  same  assump¬ 
tions  on  F, 

■>x  k 

-’0  u  F(u>  du  “  *  CTkl  F(o)  +  Tk2  F(Jkx)  +  Tk3  F(x)3  , 

(2.4-15) 
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where  the  expressions  for  the  barred-quantities  in  (2.4-15)  are 
determined  from  (2.4-13)  by  substituting  for  the  expressions 


(2.4-16) 


The  evaluation  of  the  defined  in  (2.4-14)  can  be  performed 
approximately  by  an  integration  by  parts  to  obtain 


wk  =  CTln“M 


Now  the  function  ^  ~  is  an  analytic  function  of  u,  and  the 
assumption  that  it  is  slowly  varying  allows  application  of 
(2.4-15).  Thus: 


(2.4-17b) 

By  breaking  up  the  interval  of  integration  in  (2.4-10)  into  two 
parts:  a  <  u  <  u  and  u  <  u  <  b  and  applying  (2.4-11)  - 
(2.4-17),  the  values  of  the  moments  are  determined,  subject  to 
the  assumption  that  the  R.  and  the  S.  vary  no  faster  than  cubics 
over  a  half  mesh  point  spacing. 
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3.  ANALYSIS  OF  AN  ELLIPTIC  DOME  ABOUT 
A  CIRCULAR  TRANSDUCER 

The  dome- transducer  model  described  in  this  section 
is  shown  in  Fig.  3-1.  The  transducer  is  described  by  the  circle 
r  =  a  in  polar  coordinates  r  and  9 .  The  dome  lies  along  the 
curve  §  =  §o  in  elliptic  coordinates  z  and  rj .  The  pressure 
field  between  the  transducer  and  dome  will  be  denoted  by  p^ 
and  that  exterior  to  the  dome  by  P2* 

At  the  transducer  surface  a  fluid  particle  velocity 
distribution  is  specified  to  simulate  the  velocity  distribution 
across  the  active  array.  This  velocity  distribution  is  taken 
to  be  two-dimensional,  that  is,  the  velocity  v^  of  the  elements 
of  the  jth  stave  is  assumed  to  extend  indefinitely  along  the 
axial  direction  of  the  circular  cylinder.  The  mathematical  rep¬ 
resentation  of  this  velocity  distribution  is 


v 


(3) 


lO  ,  elsewhere. 


(3-1) 


Here  9^  is  the  angular  location  of  the  center  of  the  jth  stave 
and  a  is  half  the  stave  angular  width.  This  distribution  may  be 
expanded  in  a  Fourier  series  [7j, 


Vrav  =  e  ^  6^  cos  n9  +  6s  sin  n9  .  (3-2) 

)  —  n  n 

n=0 


The  solution  p^  of  the  wave  equation  may  be  represented 
in  several  ways.  In  order  to  match  the  boundary  condition  at 
r  =  a  it  is  convenient  to  use  the  Bessel  function  expansion, 
which  is  characteristic  of  circular  cylindrical  coordinates: 
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CO 

P1  =  L  [An  Jn(kr)  +  Bn  Nn(kr)]  cos  ^ 
n=0 

+  Jn(kr)  +  Nn(kr)l  sin  n3 


(3-3) 


Here  J  and  N  are,  respectively,  the  Bessel  and  Neumann  functions 

^  ^  CSC 

of  order  n  and  argument  kr.  The  unknown  coefficients  A  ,  A  ,  B 

°  n  n  n 

and  Bs  are  to  be  determined  by  matching  boundary  conditions  at 
the  transducer  and  at  the  dome.  The  radial  component  of  the 
fluid  velocity  is  related  to  the  pressure  field  by 


This  relation, 


and 


1  3  p 

v  =  -i —  r-*-  . 
lujp  3  r 


in  conjunction  with  Eqs.  (3-2)  and  (3-3),  implies 


An  Jn(ka)  +  n  Nn(ka>  =  iPc  6 


n 


(3-4) 


Aj  Jn(ka)  +  B*  Nn(ka)  =  ipc  6 


n 


(3-5) 


The  primes  on  the  Bessel  and  Neumann  functions  designate  dif¬ 
ferentiation  with  respect  to  argument. 

According  to  Assumptions  7  and  8  in  Section  1,  a  pres¬ 
sure  difference  &p(rj)e  1U)t  across  the  surface  produces  a  normal 
vibrational  displacement  w(^)e  :Lujt  which  is  related  to  the  pres¬ 
sure  difference  by 
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'-pH 


w(n)  . 


(3-6) 


Here  ag  is  the  shell  surface  mass  density. 

General  solutions  to  the  wave  equation  may  be  written 
for  the  pressure  fields  inside  and  outside  the  shell  in  terms  of 
elliptic  coordinates: 

CO 

pl  ’  1  [aSCen<!>  +  Bn  Fe7n<?>l  C8„W 
n=0  J 

+  [an  Sen('!)  +  3n  Geyn(l)j  sen(ri)  ’  (3~?) 

co 

p2  “  'n  Men<5)  «„(l)  +  v°  Nen(S)  se  (,)  .  (3.8) 

n=0 

In  the  above  equations,  cen(rj)  and  sen(7?)  are  Mathieu  functions 

and  Ce  (?)»  Feyn(-)>  Sen(?)j  and  Geyn(?)  are  associated  Mathieu 

functions.  The  Hankel  type  associated  Mathieu  functions  in  Eq. 

(3-8)  are  defined  by  Me  =  Ce  +  i  Fey  ,  Ne  =  Se  +  i  Fey  . 
v  J  n  n  Jn’  n  n  7n 

These  functions  insure  that  there  i.j  only  outgoing  radiation  in 

6  0  6  0  6  O 

regions  outside  the  dome.  The  a  ,  a  ,  3  ,  B  ,  Y  and  v  are  un- 
known  coefficients  which  are  to  be  determined  by  matching  boundary 
conditions  at  the  dome  and  transducer. 

At  the  inner  and  outer  surfaces  of  the  dome  the  fluid 
velocity  component  normal  to  the  dome  is  equal  to  Matching 

these  boundary  conditions  at  the  dome  middle  surface,  q  =  5  , 
gives 


1 

=  oW  _  1 

(*P2\ 

ixo 

Un  .W 

*  t  imo 

Un  J 

-  o 


(3-9) 
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In  elliptic  coordinates  the  normal  derivative  is  computed  accord¬ 
ing  to 


ti„  ■ 


o  h  V  sinh  +  sin2ri 


M„  ■ 


(3-10) 


where  h  is  half  the  interfocal  separation.  Combining  Eqs.  (3-6), 
(3-9),  and  (3-10)  yields 


EPl  -  P2J___  =  “ 


physinh'T  +  sin  n 
*  o 


sr&L- 


(3-11) 


i  -ipn  Cen(5e>  +  6n  Feyn(?o> 
n=0  L 


-  Yn  Men(?o>]  cen(,,) 


+  an  Sen(5o>  +  Sn  Gey„(?o>  *  ~<n  Nen<?o>]  sen(,4 

03  — ' 

=  - j  3  ^  L  Y^4(Vs%(^>4<?o)sen(T|)- 

phV  sinh^  +  sinzr|  n=0 

°  (3-12) 

where  the  primes  on  the  Men  and  Nen  designate  differentiation  with 
respect  to  argument.  There  will  be  coupling  among  the  ce^ (y)  and 
sen(V)  modes  unless  og  is  made  to  vary  asVsinh  §  +  sin  77.  This 
variation  in  <j  can  be  considered  as  resulting  from  a  dome  thick- 

O 

ness  which  varies  as  the  distance  between  confocal  ellipses. 

Assuming  such  a  variation  of  o  and  employing  the  linear  indepen- 

s 

dence  of  the  Mathieu  functions  leads  to  the  relations 
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an  Cen<?o>  +  FeV'o>  '  Y„  Men(S0)  =  -  ^  v'  Men(?o)  ,  (3-13) 
an  Sen<^-o)  +  Sn  Geyn(-o^  '  Yn  =  ‘  j3T  Yn  Nen(?o^  ■  (3-14) 


o  .  e 


I  2  2~ 

where  a  =  a  Vsinh  5  +  sin  ij.  The  condition  in  Eq.  (3-9)  may 

s  o 

be  simplified  by  using  Eqs.  (3-7)  and  (3-8)  and  the  linear  in- 


n  n'"oy 


n  n'-o' 


'S 


n'-  o' 


n  •'n'  o' 


Thus , 

Yn  MeX  • 

(3-15) 

Yn  <«o>  • 

(3-16) 

Equations  (3-13)-(3-15)  may  be  used  to  express  ae,  a°. 


0^,  and  0n  in  terms  of  y^  and  y°  . 

an  =  f1  -  §  Me(S  ) 
n  [_  w  nv  oy 

Sn=  |_i+CT«e>„> 

an  i1  -  §  NeX> 

3°  =  [i  +  §  Ne#  (?  ) 
n  |_  w  nv  o' 


These  results  are 

F%(?o>]  Yn  .  > 

(3-17) 

CeX]  Yn  ’ 

(3-18) 

GeyX>]Y°  ’ 

(3-19) 

SeX>l  Yn  • 

(3-20) 

Here  W  is  the  Wronskian  of  the  associated  Mathieu  functions  and 

R  =  °7-  .  For  the  definition  of  these  functions  used  in  the  final 
oh 

computations,  W  =  2/ti. 

At  this  point,  the  pressure  field  p^  inside  the  dome 
has  been  represented  in  both  circular  and  elliptic  coordinates. 
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The  two  representations  must  be  identically  equal  at  all  points 
within  the  dome.  A  sec  of  equations  for  determining  the  unknown 
coefficients  y®  and  can  be  found  by  demanding  that  the  two 
representations  of  p-^  be  equal.  This  procedure  is  facilitated 
by  employing  the  expansions  of  the  elliptic  coordinate  eigen¬ 
functions  in  terms  of  the  circular  coordinate  eigenfunctions. 
These  are: 


Cen(?)  cen(q) 


'L  (~l)m+q  Ap  Jp(kr)  cos  p 
m=0 


p  =  2m,  n  even 
p  =  2mfl,  n  odd  (3_2i) 
q  =  n/2,  n  even 

[q  =  ,  n  odd 


Fey  n  (5)  cen(q)  =  (-l)nH'q  Ap  Np(kr)  cos  p  0  , 

m=0 

oc 

Sen(§)  s  en (h )  =  ^  (-l)mfq  Bp  Jp(kr)  sin  p  9  , 

m=0 
00 

Geyn(?)  sen(ri)  =  L  (-1)^  Bp  Np(kr)  sin  p  6  . 

m=0 


(3-22) 

(3-23) 


(3-24) 


The  An  and  B*}  are  the  Fourier  coefficients  in  the  expansions 
P  P 


ce^Cl)  =  2,  Ap  cos  P  ^  » 


m=0 


sen(ri)  = 


00 

=  V  Bn 


m=0 


sin  p  p  . 
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By  equating  the  expressions  in  Eqs.  (3-3)  and  (3-7) 

and  substituting  Eqs.  (3-21) - (3-24)  for  the  elliptic  coordinate 

CCS  s 

eigenfunctions,  the  coefficients  An,  Bn,  An,  and  Bn  are  found 
in  terms  of  a^,  0^,  a°,  and  0°.  Then  Eqs.  (3-4),  (3-5),  (3-17), 
(3-18),  (3-19),  and  (3-20)  may  be  used  to  eliminate  all  the  coeffi¬ 
cients  except  the  Ya  and  y°  .  The  final  equations  are  shown  in 
Table  3-1. 

The  equations  in  Table  3-1  constitute  four  infinite 

p  o 

sets  of  linear  equations  in  the  y  and  Yn*  Numerical  experimen¬ 
tation  indicates  that  they  are  convergent  in  the  following  sense. 
Suppose  a  set  of  2M  of  the  yn  and  2M  of  the  yn  are  obtained  by 
solving  the  first  M  equations  of  each  set,  truncating  the  sum  in 

p 

each  equation  after  the  M  term.  If  the  first  N  of  the  Y2n  , 

Y2n+1*  Y2n+1’  ancl  Y2n+2  not  c^an8e  in>  saY>  first  six 

significant  figures  as  M  is  increased  past  some  value,  denoted 

by  Mn,  then  the  equations  are  said  to  have  converged  to  N  terms 

for  an  by  set  of  equations.  When  a  sufficient  number  of 

ye  and  y°  have  been  computed  in  the  above  manner  for  the  series 
n  n 

in  Eq.  (3-8)  to  converge,  the  problem  solution  is  completed. 
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4.  THREE-DIMENSIONAL  CYLINDER-DOME  MODEL 


The  mathematical  model  described  in  this  section  (Fig. 
4-1)  is  an  extension  of  the  two-dimensional  cylindrical  trans¬ 
ducer  and  concentric  shell  model.  In  the  present  model  the  active 
region  of  the  transducer  is  finite  in  height.  This  generaliza¬ 
tion  of  the  model  permits  the  computation  of  three-dimensional 
beam  patterns  and  directivity  indices  of  cylindrical  arrays. 

The  pressure  field  of  a  singl..  rectangular  element  will 
first  be  obtained.  Then  superposition  will  be  employed  to  deter¬ 
mine  the  array  beam  pattern.  Let  the  particular  element  chosen 
be  such  that  the  x-z  plane  passes  through  its  center.  The  time 
dependent  factor  (e  lat)  will  be  suppressed  in  the  following. 


The  velocity  distribution  on  the  cylinder.  r=a,  is  non¬ 
vanishing  only  on  the  element  surface.  That  is. 


v  = 


|  v0  |z  -  zc\  <  ZD,  M  < 
1 

l 0  elsewhere 


(4-1) 


where  z  is  the  z  coordinate  of  the  element  center,  z  is  the 
half  height  of  the  element,  and  is  the  half  angular  width  of 
the  element.  It  is  convenient  to  take  the  Fourier  transform  of 
all  the  fields  which  depend  on  the  z  coordinate. 

The  Fourier  transform  of  v  is  defined  to  be 


!jl=  Lv  e"uz  dz  m<*( 


V  = 


. 


1 

M  >  J 


(4-2) 
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FIG.  4-1 -GEOMETRY  FOR  A  THREE  DIMENSION  CYLINDRICAL 

TRANSDUCER  RADIATING  INSIDE  A  CONCENTRIC  DOME 
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The  Fliigge  shell  equations  are  assumed  tc  describe  the  be¬ 
havior  of  the  dome  [81.  The  following  dynamic  equations  are  Fliigge’ s 
equations  with  an  added  D’Alembert  force.  The  fluid  is  inviscid 
and  therefore  there  are  no  external  forces  acting  in  the  s  or  z 
directions : 


(4-11) 


The  quantities  in  Eqs.  (4-9)  -  (4-11)  are:  b,  the  middle  radius 


of  the  dome; 
material;  h, 


v,  Poisson's  ratio;  p  ,  the  density  of  the  dome 

s  *  ti¬ 
the  thickness  of  the  dome;  k  =  - D  =  E  h;and  u  ,u  , 

12bz  s  z  ® 
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and  ur,  the  displacement  components  of  the  middle  surface  of  the 
dome  in  the  z,  r,  and  r  directions,  respectively. 

Equations  (4-9)  -  (4-11)  must  be  Fourier- transformed.  It 
is  assumed  that  the  displacements  may  be  expressed  as  Fourier 
series. 

OD 

uz=  X  cos  03  *  (4-12) 

n=o 

OD 

u^  =  V  vn(z)  sin  rco  ,  (4-13) 

n=o 

OD 

ur=  2  Wn(z)  cos  r£p  *  (4-14) 

n=o 


The  Fourier  transforms  of  u  u  ,  and  u  are  defined  respectively 

z  cp  r 

as 


u 


z 


u 


z 


(4-15) 


Ucp 


1 


J 


u 


-iXz 


dz 


CP 


(4-16) 


u  =  ~r  1  u  e  dz  . 

r  \/2tt  r 


(4-17) 


The  displacement  components  of  the  dome,  Eqs.  (4-12)  -  (4-14),  can 
now  be  Fourier- transformed : 
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u,  =  ^  Un  (X)  cos 


(4-18) 


^  Vn  (X)  sin  n c? 


(4-19) 


5r  =  I  «n  (X> 


cos  rep  . 


(4-20) 


It  is  now  convenient  to  Fourier- transform  Fliigge's 
shell  equations.  In  addition  to  Fourier- transforming  Eqs.  (4-9)  - 
(4-11),  the  transformed  Fourier  series  representation  for  the  dis¬ 
placements,  Eqs.  (4-18)  -  (4-20),  are  substituted  into  Fliigge's 
equations.  To  reiterate  Assumption  4,  the  displacement  components 
have  harmonic  time  dependence.  The  relationships  between  trans¬ 
formed  displacement  components  and  transformed  velocity  com¬ 
ponents  of  the  shell  are  therefore: 


v„  -  -it  u  , 
z  r  ’ 


V  =  -IX’  u  , 
cp  * 


V  =  -IT  u 
r  r  , 


where  vz,  v  ,  and  vf  are  respectively  the  transformed  velocity 
components  in  the  z,  co,  and  r  coordinate  directions.  The  trans¬ 
formed  shell  equations,  in  terms  of  velocity  components,  can  be 
written  as  follows: 
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\  ["x2b2  "  n2  (1+k*) 


On  [-i>bn(i^]  +  Vn  [-n2-X 
wn  [-n  -k*X2 


1 2b2o  h  - 
s 

]  +  Vn  £iXbn' 

D 

*(ib3X3  - 

V 

iXbn2  <*?>) 

2b20f> 

-  k*X 2  j  b2 

■b2n^] 

=  0  V 

'»<"»]  + 

Cn  [I'  Hn(Kb 

(1+v) 


')]  - 0  • 


(4-21) 


,2  2, 

b  x  Psh 

D 


(4-22) 


+ 


5»[ 


\  [n+k*X2b2n(^.)]  +Wn 


iXbv  -  k*iXbn^ 


0  Ck^)  +  k*  ix3b3 


1  +  k*  (X4b4  +  2\2b2n2  +  n4)  - 


k*  (2n2-l)  -  £  Pghx2 


=  0  . 


(4-23) 


To  obtain  Eq.  (4-23)  from  Eq .  (4-11),  it  was  necessary 
to  assume  solution  forms  for  the  pressures  px  and  p2.  The  pressure 
in  region  1  (see  Fig.  4-1)  is  denoted  by  p^  and  in  region  2  by 
p  When  Eq.  (4-11)  is  Fourier- transformed,  the  pressures  p1  and 
p^  transform  into  p-^  and  p2,  respectively.  The  transformed 
pressures  pl  and  p?  must  satisfy  the  transformed  Helmholtz 
equation,  Eq.  (4-6).  The  form  of  the  solution  to  Eq .  (4-6)  is 
known  and  can  be  written  as 
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-Y 


~AnJn(Kr)  +  BnNn^Kr^  cos  025  > 


(4-24) 


h  =  £  C^H^CKr)  cos  no 


(4-25) 


It  is  necessary  to  require  that  K  M  to  insure  that 
the  Bessel  function  of  the  second  kind  in  Eq.  (4-24)  remains 
finite  in  the  range  of  interest. 

One  boundary  condition  which  must  be  satisfied  is  the 
transducer-fluid  interface  condition  at  r=a,  Eq.  (4-8).  The  pres¬ 
sure  is  given  by  p^,  and  v  is  given  by  Eq.  (4-3).  The  transformed 
boundary  condition  at  r*a  yields  the  relation, 


v;<Ka> +  wKa>  * 


i'J>Pl/2  c  v  zcp„  sin  Xz„  sin  no. 


n  o  o’o 


K(n) 


>  (4-26) 


where  the  prime  denotes  the  differentiation  with  respect  to  argu¬ 
ment,  and  e  is  the  Neumann  symbol  e  =  1,  e  =  2  n  >  o.  Two 
more  boundary  conditions  must  be  satisfied  at  the  middle  surface 
of  the  shell.  The  radial  velocity  of  the  middle  surface  of  the 
shell  must  equal  the  normal  gradient  of  the  pressure  evaluated 
at  r=b.  This  requirement  applies  to  representations  of  the 
pressure  both  inside  and  outside  the  shell. 

Satisfaction  of  the  two  shell  middle  surface  boundary 
conditions  requires  that 


uu2p  W 


Vn(Kb)  +  BnNn(Kb)  =  iC 


(4-27) 


(Kb)  =  ac  P  Wr 


(4-28) 
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The  solution  is  now  such  that  there  are  six  equations 

and  six  unknowns.  Equations  (4-21)  -  (4-23),  (4-26),  and  (4-27)  - 

(4-28)  are  the  six  equations,  and  A  ,  B  ,  C  ,  U  ,  V  ,  and  W  are 
v  n  ’  n*  n  n  n  n  n 

the  unknowns.  After  tedious  algebra,  the  unknown  coefficients 
Cn  may  be  obtained: 


cn(x) 


n 


m.  hE  ttkV1)  (Kb)  A 

rW  (Ka)  +  — i 


n 

TT- T“ 

2ba’  o 


- 5 - ? 

n  l  2m2ra6m3  -  m3m5  + 

m9  +  - ’5 - 


(m^in^  i  m2) 


(4-29) 


where 

ia’p  'Jl  £  v  z  co  sin  Xz^  sin  n  co 
_  '  n  o  o  o  _ o  _ o 

‘  "  =  K  (*)in  Xzo  n  %  ’ 

\  -  j;(Ka)  N^(Kb)  -  j;(Kb)  N’ (Ka)  , 

2  2 

o  o  o/i  \  (i)  p„b  h 

n»1  — X2b2  -  n2  (1  +  k +  - D -  ’ 


m 


2 


iXbn 


) 


m3  =  iXbv  +  k*  [iX^b^  -  iXbn2^ 


(V)]  > 


.2  ,  2,  2  /1-v' 


2  3  u2, 


b2o)2p  h 
s 


m3  =  -  n‘  -  X  ‘‘b4’  ^ — 2~j  ~  k*  ^  ^  b^(l-v)  +  - ^ 
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m 


6 


n  -  k* 


2  2 
X  d  n 


> 


IHr 


=  1  +  k*  '*  2l2  +  n2)2  -  k*  -  1) 


b2o  hu)2 
s 

- B - 


With  the  coefficient  C  thus  obtained,  the  transformed 

n 

pressure  in  region  2  is  determined.  The  problem  to  this  stage 

has  been  even  about  the  x-z  plane.  If  cd  f  0  as  assumed,  but  if 

instead  the  element  is  located  at  cd  then  cd  -*(cd  -  cd  )  : 

c  c 

cos  nm  -  cos  n(cp  -  cp£) 

cos  n(cp  -<?>)  =  cos  np  cos  ncp  4-  sin  np  sin  no  . 
c  c  c 


P2 


Equation  (4-25)  can  be  rewritten  as, 
00 

id) 


[CnH^L‘/  (Kr)cos  ncpc  cos  np  +  cnH^2^(Kr) sin  ncpc  sin  np]#  (4-30) 


n=o 


The  inverse  Fourier  transform  of  P2  yields  a  formal 
solution  for  p: 


00 

Cn(X)[ccs  ncpc  cos  ncp  +  sin  ncpc 

n=o 


sin  ncp]  e  IP  '  (Kr)  dX 


(4-31) 


The  pressure  p,  as  written,  is  in  terms  of  cylindrical  polar  coor¬ 
dinates.  It  will  be  convenient  to  express  p  in  spherical  coor¬ 
dinates  (R,6,cd).  From  Fig.  (4-1),  r  =  R  cos0  and  z  =  R  sin0. 


P 


»  CO 


,  _  ^XRsinO 

C  (X)[cosnp  cosncp  +  sinnp  sinnp]e  H 

n.  c  c 

.00  ^ 
n=0 


(i) 

n 


(KRcos0 ) dX. 
(4-32) 


73 


r 

ft 


6500  TRACOR  LANE  A’jST'N  TEXAS  78721 


The  crux  of  the  problem  is  the  evaluation  of  the  inte¬ 
gral  in  Eq.  (4-32).  The  method  of  steepest  descent  can  be  applied 
to  Eq.  (4-32)  [9].  The  method  of  steepest  descent  is  a  technique 
to  tind  an  asymptotic  approximation  to  a  certain  class  of  inte¬ 
grals,  namely,  integrals  of  the  form 

I(r)  =  J  g(z)  er£^dz  , 
c 

where  r  is  large  and  positive.  By  applying  steepest  descent, 
it  has  been  tacitly  assumed  that  only  farfield  pressure  results 
will  be  forthcoming.  Equation  (4-32)  must  be  of  such  a  form  that 
the  contour  of  integration  can  be  deformed  so  that  an  asymptotic 
expansion  of  the  Hankel  function  can  be  made.  Further,  the  con¬ 
tour  must  be  deformed  to  pass  through  the  saddle  points  of  f(z). 
Equation  (4-32)  meets  all  necessary  requirements  for  a  steepest 
descents  analysis.  The  results  are, 


p(R»  9  ,  cp)  =  2- 


ikR 


R 


I 


n=o 


-i-^(nfl) 


Gn(s>  V 


(4-33) 


where 


Gn(co,  =  Gn^o^cos  ^c  cos  ^  +  s^n  ^c  s:''n  >  anc* 

\j2jr 

\  =  k  sin  cp  . 

o 

The  farfield  pressure  is  now  completely  determined. 


6500  TfiACOR  LANE,  AUSTIN.  TEXAS  78721 


5.  3PHERICAL  TRANSDUCER-DOME  MODEL 

The  analysis  of  the  spherical  transducer  and  concen¬ 
tric  spherical  dome  was  presented  in  a  previous  technical  memoran¬ 
dum  [10] .  The  earlier  analysis  will  be  outlined  in  this  section, 
then  the  extension  of  the  analysis  necessary  to  compute  mutual 
acoustic  impedances  will  be  presented. 

The  geometry  of  the  spherical  transducer  and  spherical 
dome  is  shown  in  Fig.  5-1.  The  first  step  in  the  analysis  is  the 
solution  for  the  sound  pressure  radiated  by  a  single  piston.  The 
vibrating  piston  will  produce  a  sound  pressure  field  between  the 
dome  and  transducer.  This  pressure  field  is  a  solution  to  the 
wave  equation  and  can  be  represented  by  the  expansion, 


P1  =  lAmin>(kr>  +  Vm<kr>]  Vcos  ♦)  ■  (5-1) 

m=0 

This  form  is  standard  for  solutions  of  the  wave  equation  in 
spherical  coordinates  (r,  )  if  there  is  symmetry  about  the  axis 

from  which  the  angle  i'  is  measured.  The  Pm  are  Legendre  poly¬ 
nomials  of  argument  cos  and  the  j  and  n  are,  respectively, 

the  spherical  Bessel  and  Neumann  functions  of  argument  kr.  The 
A  and  B  are  constants  which  can  be  adjusted  to  make  the  solu- 
tion  form  fit  whatever  boundary  conditions  are  prescribed  on  the 
transducer  and  shell  surfaces.  The  boundary  conditions  arise  in 
equating  the  fluid  and  solid  velocity  components  normal  to  the 
fluid-solid  interfaces.  Velocity  components  tangential  to  the 
interfaces  are  not  necessarily  equal  since  the  fluid  is  inviscid. 

For  the  case  at  hand,  the  fluid  velocity  is  equal  to 
the  piston  velocity  on  the  piston  surface  and  is  zero  on  the 
remaining  transducer  surface.  The  piston  velocity,  v  e"lat, 
will  be  considered  to  be  purely  radial,  so  the  boundary  condition 
at  the  transducer  surface  r  =  a, is 
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00 


m=0 


/ 


(ka)  +  B  n'  (ka) 
'  m  mv  ' 


P 


m 


(cos  v) 


JVo,  0  *  ,  <  Q 
1 0,  a  <  t  *  ~ 


(5-2) 


Here  the  angular  extent,  of  the  piston  is  0  <  v  <  a,  the  radial 
fluid  velocity  component  being  zero  outside  of  this  range.  The 
primes  on  the  j  and  n  denote  differentiation  with  respect  to 
argument . 

Whi?_e  the  pressure  field  between  the  transducer  and 
dome  has  the  form  given  in  Eq.  (5-2),  the  sound  pressure  field  p2 
outside  the  dome  must  have  a  form  which  describes  only  outgoing 
radiation  for  large  r.  This  form  is 


P2  - 


I 

m=0 


m 


hm(kr) 


P  (cos  *) 
m 


(5-3) 


where  h  =  j  +in  is  a  spherical  Hankel  function  and  the  C 
are  determined  by  the  boundary  condition  at  the  outer  dome  sur¬ 
face,  Since  the  dome  is  thin  relative  to  an  acoustic  wavelength, 
it  is  sufficient  to  apply  the  velocity  continuity  condition  at 
the  dor  i  midsurface,  r  =  b,  rather  than  at  the  actual  fluid- 
dome  interfaces.  Thus, 


1 

i'xc 


_l_lP2 

iwc  5  r _ 

r=b 


(5-4) 


w  being  the  radial  velocity  distribution  of  the  dome.  The 
Legendre  polynomials  form  a  complete  set  over  the  range 
0  <  ;  ",  hence  the  dome  velocity  admits  of  an  expansion 
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✓ 


w  =  Y  Mf  P  (cos  t)  .  (5-5) 

/  .  mm7 

m=0 

Equations  (5-2)  and  (5-5)  give  three  conditions  on  the  four  sets 

of  unknowns,  A  ,  B  ,  C  and  W  .  The  fourth  condition  must  corne 
m’  m  m  m 

from  an  analysis  of  the  motion  of  the  dome. 

In  the  earlier  report  [10],  it  is  shown  that  the  appli- 

—  T  Jl*  t" 

cation  of  a  load  distribution  of  the  form  a  P  (cos  £)  e  to 

mm 

the  dome  will  cause  a  velocity  distribution  w  -  ^  p  (cos  i)  e~^aC 

£  m 
m 

where  is  the  modal  impedance  of  the  dome.  The  driving  force 
on  the  dome  in  the  problem  under  consideration  is  just 


Pi  ~  P2 


■I 


AJm(kb) 


Bn  (kb) 
mm 


C  h  (kb)  , 
mm  j 


Pm(cos  *) 


m=0 


(5-6) 


Thus  the  dome  response  is 


w=Y  ~  A  j  (kb)  +  B  n  (kb) 
L. i  Z  mJm  '  m  mN  7 


m=0 


m 


C_h^(kb)j  Pm(cos  £•) 


m  m 


m 


(5-7) 


Equation  (5-7)  is  the  fourth  condition  required  to  solve  for  the 


m : 


B 


m 1 


m: 


and  W 


m 


The  solution  for  the  expansion  coefficients  is  straight 
forward,  involving  only  the  use  of  the  orthogonality  of  the 
Legendre  functions  to  derive  simple  algebraic  equations  from 
which  the  expansion  coefficients  can  be  obtained.  For  farfield 
radiation  characteristics,  only  the  coefficient  C  is  required. 

It  is  given  by 
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P  t (cos  a)  -  P  , , (cos  a) 
m-1  '  m-f-1  ' 

Z  “ 

h'(ka)  +  -r^— (kb)  h'  (kb)r  j  '  (ka)  n  (kb)'  -  j  '  (kb)  n'(ka)'' 
mx  '  iccv  '  mv  1  mv  J  Jmv  '  mv 

(5-8) 


To  compute  the  farfield  distribution,  the  asymptotic 


form  of  hm(kr)  for  large  r. 


,  •  /,  ,  l'"1'1  ikr 

Lxm  hm(kr)  =  e 

kr  -  ® 


can  be  substituted  into  Ec.  (5-3).  Thus, 


Lim  p(r,c)  =  pcv 


i(kr-a’t) 


F(cos  •£)  , 


(5-9) 


r  -  “ 


where 


F(cos  •))  = 


Ik  ) 


i"m  Pm_1(cos  a)  "  Pm+1  (cos  a) 


m=0  h' (ka)  +  -» — (kb)  h'(kb)[j  (ka)n' (kb)-j' (kb)n' (ka)l 
mN  '  10  c'  '  m  /  Jm  7  m  7 


This  expression  determines  the  farfield  characteristics  for  a 
single  piston. 

The  farfield  sound  pressure,  p(r,9,§),  of  an  array  of 
pistons  can  be  found  by  superposition  of  the  fields  produced  by 
the  individual  pistons.  The  piston  whose  center  is  at  9j,  on 

the  transducer  has  velocity  Vje~lat:.  The  cosine  of  the  angle 
between  8  ,  §  and  9j,  is  just 
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cos  *j  =  cos  3  cos  3j  +  sin  3  sin  9  T  cos  (§  -  ? j)  .  (5-10) 

Superposition  of  the  pressure  fields  of  an  array  of 
N  pistons  gives 


i(kr  -  cut) 

P(r,0  ,$)  =oc^ - - -  D(0,*)  ,  (5-11) 


where 


N 

D(0  ,$)  =  ^  vj  F(cos  $j)  . 
J=1 


The  farfield  intensity  distribution  I  is  p*  p,  where  the 
asterisk  denotes  the  complex  conjugate.  Therefore, 


I  =  ££  -ijr  |  D(0  ,$  )  |  ^  .  (5-12) 

r 

The  pressure  field  p-^  for  a  single  piston,  evaluated 
at  the  transducer  surface,  is  found  to  be 

P1  =  Vo  ^^rn  [Pn-l(cos  a>  "  Pn+l(cos  Pn(cos  *  )  » 

n  (5-13) 

where 

K  9 

hn(ka)  +  ^~(kb)  Zhn(kb)[  jn(ka)n^(kb)  -  j^(kb)nn(ka)j 
n  h^(ka)  +£^2.(kb)  2h^(kb)(  j^(ka)  n^(kb)  -  j^(kb)n^(ka)] 

This  pressure  fieid  causes  a  net  force  on  the  other  elements  in 
an  array. 
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By  use  of  the  addition  theorem  for  Legendre  polynomials 
this  solution  may  be  represented  in  a  different  polar  coordinate 
system,  whose  z-axis  passes  through  the  center  of  a  second  piston. 
Then  the  total  force  on  the  second  piston  due  to  the  pressure 
field  of  the  first  may  be  determined  by  integrating  over  the  sur¬ 
face  of  the  second.  The  force  on  the  second  piston  per  unit 
velocity  of  the  first  is  called  the  mutual  acoustic  impedance. 

If  the  elements  have  an  angular  separation  9 ,  the  expression  for 
the  mutual  impedance  is 


9  "  [P  ,(cos  a)  -P  ,(cosa)]2 

=  ina  oc)  T  — - — -  P  (cos  S)  . 

/  n  « '  / 


2n+l 


(5-14) 


For  5  =  0,  this  gives  the  self  impedance  of  the  piston. 


nt  uwiwim  ..hwpwuw  iw'wwii  <*■•»» 
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6.  THICK  DOME  MODEL 

The  mathematical  model  illustrated  in  Fig.  6-1  consists 
of  an  infinite  length  cylinder  and  concentric  dome.  A  velocity 
distribution  independent  of  the  axial  direction  is  specified  on 
the  active  surface  of  the  transducer.  It  is  desired  to  fj.nd  the 
resulting  sound  pressure  field  in  the  fluid  in  regions  1  and  3 
of  Fig.  6-1.  The  boundary  value  problem  is  a  three  region  prob¬ 
lem:  two  fluid  regions  and  the  elastic  dome  region. 

The  general  solution  to  the  scalar  wave  equation  in 
region  1  is 


P 


1 


33 

r  1  0 
=  *-AnHn^kr^  +  B  H;(kr)]cos  ne  , 
n=0  n  n 


(6-1) 


where  An  and  Bn  are  constants  to  be  determined  by  boundary  condi¬ 
tions.  Equation  (6-1)  involves  only  cosine  terms  as  the  problem 
to  be  analyzed  will  arbitrarily  be  made  even  about  0=0°. 

In  region  3  the  second  pressure  field  is 


P 


3 


\  1 

:  _Q  CnHn(kr)C0S  n9 * 


(6-2) 


?  Equation  (6-2)  satisfies  the  Sommerfeld  radiation  condition, and 

,  Cn  is  a  constant  to  be  specified  by  the  boundary  conditions. 

r  At  the  transducer-fluid  interface,  the  radial  component 

!  of  the  fluid  velocity  must  equal  the  velocity  of  the  transducer. 

|  The  radial  velocity  of  the  transducer  may  be  expanded  in  a  Fourier 

|  series  of  the  form. 


CO 

v  -4-  V  6 
r  ioc  _n  n 
n-U 


cos  no 


(6-3) 
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The  boundary  condition  at  the  transducer-fluid  interface  may  now 
be  matched  to  yield, 


A 

n  n 


(ka)  +  B  H^2) ' 
n  n 


(ka) 


(6-4) 


where  the  prime  (' )  denotes  differentiation  with  respect  to  the 
argument  of  the  Hankel  function. 

A  means  must  now  be  devised  to  describe  the  motion  of 
the  thick  shell.  Let  v  be  the  velocity  field  vector  of  the  thick 
shell,  and  u  be  the  displacement  field  vector  of  the  shell.  As¬ 
sume  u  and  v  are  harmonically  time  dependent.  The  velocity  and 
displacement  field  vectors  are  related  as  follows: 


v  =  -ia'u  .  (6-5) 

An  important  property  of  a  vector  is  that  it  tray  be  represented 
as  the  gradient  of  a  scalar  plus  the  curl  of  a  vector [11].  The 
vector  v  may  be  written  as 


v  =  W  +  V  x  if  . 


(6-6) 


For  this  two-dimensional  problem,  $  =  tlr k .  Both  cp  and  il<  satisfy 
the  scalar  Helmholtz  equation,  while  v  is  a  solution  of  the 
vector  Helmholtz  equationC  12]  .  The  assumed  solution  forms  for 
cp  and  i!i  must  satisfy  the  Helmholtz  equation;  and  because  the 
problem  under  consideration  was  assumed  even,  the  solution  forms 
must  yield  a  radial  velocity  component  which  is  an  even  function 
of  Q  and  a  tangential  velocity  component  which  is  an  odd  function 
of  8.  If  the  following  forms  are  assumed  for  cp  and  if, 


LANE.  AUSTIN.  TEXAS 


78721 


T  =  I  (c  r)  -f  EnH^2'  (r  r)]  cos  n  S  ,  (6-7) 

n=0 

CO 

*  =  2:  [ F  Hp  (§ r)  +  G  H^2)(§r)]  sin  n  9  ,  (6-8) 

n=0  n  n 


then  co  and  i  satisfy  the  scalar  Helmholtz  equation  and  meet  the 
above  requirements  on  the  radial  and  tangential  velocity  com¬ 
ponents.  The  unknown  coefficients  D  ,  E  ,  F  ,  and  G  will  be 
determined  so  that  the  boundary  conditions  at  both  shell-fluid 
interfaces  are  satisfied.  The  quantities,  Q  =  w/c’  and  §  =  jj/c1', 
are  respectively  the  x^ave  numbers  associated  with  the  compres- 
sional  and  shear  waves  in  the  dome .  The  compressional  wave 
velocity  is  given  by  c 1  ,  and  the  shear  wave  velocity 

is  given  by  c ' '  •  In  the  velocity  expressions,  \  and  p 

are  the  Lame  constants,  and  0  is  the  density  of  the  shell  [13]. 

o 

From  Eqs.  (6-5)  and  (6-6)  one  finds  that 


where 


u  =  jjj-  [tfp  +  v  xf]  , 

-*  =-  -  1  3 

Sr9r  e0rS9  * 


(6-9) 

(6-10) 


The  unit  vector  in  the  radial  direction  is  and  in  the  tangen¬ 
tial  direction  e  .  If  Eqs.  (6-7) -(6-10)  are  combined,  the  re- 

o 

suits  are 


n=U 

(6-11) 


r<2) 


(I) 


,(2) 
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ur  -  i  ■  :,D ^^'CD+E  H<2)'(;r);4fiFnHW(-;r)4<;nH(2)(;r):3cosn?, 

n=0  (6-12) 

where  is  the  displacement  of  the  shell  in  the  tangential  di¬ 
rection  and  is  the  displacement  of  the  shell  in  the  radial 
direction . 

The  boundary  conditions  which  must  be  satisfied  at  the 
inner  and  outer  surfaces  of  the  thick  shell  are 


■L.  -  -  ’i  • 


(6-13) 


r]  =  “  p3  > 
Jr=e 


(6-14) 


.1=0, 

'VJ  r=b 


(6-15) 


“J  r=e 


0  , 


(6-16) 


where  the  radial  stress  in  the  shell  is  denoted  by  Srr  and  the 
tangential  stress  by  S  ^  .  In  addition  to  Eqs.  (6-13)  -  (6-16) , 
there  is  another  requirement  that  the  velocity  of  the  shell  at 
the  radii  r  =  b  and  r  =  e  must  equal  the  velocity  of  the  fluid  at 
r  =  b  and  r  =  e . 

The  equations  relating  stress  and  displacement  com¬ 
ponents  are 


„  ,  \  3u„ 

Srr  ’  <X  +  2i)  qT+F^r  +  aF^  • 


(6-17) 
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and 


S 


r0 


rlaur.aue  1  n 

uLr  TT  +  dr  r  u6] 


(6-18) 


To  obtain  the  stress  relations  in  terms  of  the  unknown  coeffi¬ 
cients  Dn>  En,  etc.,  one  substitutes  Eqs.  (6-11)  and  (6-12)  into 
Eqs.  (6-17)  and  (6-18)  to  obtain  Eqs.  (6-20)  and  (6-21)  in  Table 
6-1.  It  is  now  possible  to  satisfy  the  boundary  conditions  spec- 
fied  by  Eqs.  (6-13)  and  (6-12)  to  obtain  Eqs.  (6-21)  and  (6-22) 
in  Table  6-1.  When  Eqs.  (6-15)  and  (6-16)  are  satisfied,  the 
results  are  Eqs.  (6-23)  and  (6-24)  in  Table  6-1. 

The  two  remaining  boundary  conditions  which  must  be 
met  are  the  velocity  conditions  at  the  shell-fluid  interfaces. 

If  the  velocity  boundary  conditions  at  the  shell-fluid  inter¬ 
faces  are  matched,  and  the  unknown  coefficient  B  is  eliminated 

n 

by  use  of  Eq .  (6-4),  the  results  are  Eqs.  (6-25)  and  (6-26)  in 
Table  6-1. 


The  quantity  which  appears  in  Eqs.  (6-3)  and  (6-25) 
may  be  computed  from  the  velocity  prescription  on  the  transducer 
[7].  Equations  (6-21)  -  (6-26)  fully  describe  the  problem  under 
consideration.  These  may  be  solved  for  the  unknown  coefficients 


D  ,  E  ,  F  . 
n*  n’  n 


and  G 


nJ  n*  ~n’  ~n7  _n7  ~n 

With  all  unknown  coefficients  determined,  the  pressure 
fields  in  the  fluid  are  determined  by  Eqs.  (6-1)  and  (6-2). 
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Table  6-1 

EQUATIONS  FOR  THICK  DOME  MODEL 


(6-19) 


2 

S  =  ”  Dn:(X+2i)C2Hp//(Cr)4i^H<1>'(Cr)-^-HP(Cr)]  + 

n=0 


Enr(X+2a)C2H^2)//(Cr)  +  ^  H^2)  7  (qv)  -~-H^2)(Cr)J  + 
Fn  '(?  r)  -  iH^fir)]  + 

Gn  "Z1  -Hn2)  /(?r)  ‘  r  Hn2)^r^j  cos  n9j  * 


(6-20) 


S..-  =  — '  r  -  D  [ C hJ^1^  '  (Cr)  -  i  H^'(Cr)]  + 


n  -0  ‘ 


nr  -=  n 


iH<X> 

r  n 


E„  ~  [CHi2)/(Cr)  -  i  H^2)(Cr)3  + 


n  r  n 


Fn[?2Hn1),/(!r)  ^  (§  r)  ]  + 

Gnr  f  2H^2)  ' 7  (!  r)  -  'ip  H^2) '(“  r)  +  ^  H^z)  (?r)  ]  }  sin  n5  ) 


,(1) 


n  „(1)  (s . 


(6-21) 


Dn[(A  +  &)  C2H^1}  "(Cb)  +^H^1)/(Cb)  -^-H^>(Cb)J  + 
E/0  +  2a)  C2H<2)"(Cb)  +  VHn2),(Cb)  '  ^?Hn2)^b)]  + 
Fn%E[T  Hn1)/(rb)  "FHn1)(5b)1  + 

Gn  “E2  -  Hn2>,(;b)  '  *  Hn2)(5b)5  '  irfAnHn(1)  0*)  +  BnHn(2) 


(kb)]  . 
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Table  6-I(con't) 


(6-22) 


Dnr  (>•  +  2u)  H^fee))  + 

e 


EJO,  +  2u)  cV2)//(Ce)  +^-H^2)/(Ce)  -  ^4  H<2)  (Ce)]  + 

e 

fnT1:  ■«'(;«)  -IHW(;e)]  + 


Gn  Hn2>'«e)  -  |H<2>(|e)]  = 


1  ®  cnHn(1>(ke) 


(6‘23^iHn1)'(Cb)  '  E  HW(Cb)]^l  +  E^W'Ccb)  -  1  H<2)(Cb)]|l  + 

F^?2Ha),,(;b)  .cHa)-(;b)  +4^1)(?b)]  + 

b 

Gnl?2H<2>"(?b)  -  |H<2>'(?b)  +  4  H^2)  (§b)]  =  0 


(6-24) 


'(ie)  -  I  ^  (Ce) ]■—  +  EnljH<2>  '(Ce)  -  §  H<2>  (Ce)]^  + 

Fn[-2^1)"(?e)  -|H«'(?e)  + 

e 

Gn['2^2)"(?e)  -|H(2)'(5e)  ,,4^2)(?e))  =  0  . 
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Table  6-1  (cont'd) 


(6-25) 


Dn  1  cp:  HnX),(^b)  +  En  i  cpC  H^2)/(Cb)  +  i  cpg-  Fn  (§b)  + 


i  c-n  (2)  T  m,  H$2)/(kb)  H^'Cka) 

Gn  H ~  Hi2)(?b)  -  A  H(1>'(kb)  -  -2 - __0 - 


n  n 


7  (ka) 


Hu;  7  (kb) 

• 


(6-26) 


Dn  1  ce£  Hi1)#(Ce)  +  En  i  cpC  H^2)/(Ce)  +  Fn  — P-p-2-  H^(§e)  + 

-  Cn  Hn1}/(ke)  -  0  . 
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7.  STRUCTURE  MODEL 

The  solution  to  the  problem  of  a  cylindrical  transducer 
mounted  inside  a  concentric,  homogeneous,  cylindrical  dome  has 
been  reported  previously  [7].  The  mathematical  model  considered 
in  this  section  extends  the  solution  to  that  problem  by  incor¬ 
porating  dome  structure  into  the  mathematical  model. 

The  model  shown  in  Fig.  7-1  consists  of  an  infinite 
Cylindrical  transducer  enclosed  by  a  concentric  cylindrical  shell. 
At  locations  where  structure  elements  are  attached  to  the  shell 
there  is  added  impedance  to  the  shell  motion.  In  the  present 
analysis  this  impedance  is  assumed  to  be  mass-like,  although  the 
analysis  is  easily  extended  to  include  a  more  complicated  type 
of  impedance.  In  the  analysis  of  the  shell,  the  structure  mass 
contribution  is  included  by  inserting  an  inertial  load,  or 
D’Alembert  force,  in  the  Flugge  shell  equations. 

Flugge's  shell  equations  written  in  terms  of  displace¬ 
ments  of  the  middle  surface  are: 


3u 

375  +  ^9 


ETTT  =  0 

s 


(7-1) 


+  u 


r 


b^ 

e“K 


0  .  (7-2) 


The  radial  and  tangential  displacements  are  respectively  u^  and 

u  .  The  tangential  and  radial  force  components  acting  on  the 
9  r 

system  are  Pn  and  P  ,  respectively,  while  E  =  — ^  ,  b  is  the 
f*  r  s  i  ^ 

1  “V 

radius  of  the  middle  surface  of  the  dome,  h  is  the  thickness 
of  the  dome,  and  v  is  Poisson's  ratio.  Equations  (7-1)  and 
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FIG.  7-1  -  PLAN  FORM  GEOMETRY  FOR  A  CYLINDRICAL  TRANSDUCER 
RADIATING  IN  A  CONCENTRIC  NON-HOMOGENEOUS  SHELL 
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(7-2)  are  static  equations  but  may  be  converted  to  dynamic  equa¬ 
tions  by  including  an  additional  D'Alembert  force.  Since  the 
fluid  is  assumed  to  be  inviscid,  there  is  no  tangential  force 
acting  on  the  shell.  Equations  (7-1)  and  (7-2)  become  respec¬ 
tively  , 


J-  c  ua  + 


Esh 


9u* 


Sur 


=  0 


(7-3) 


u 


+ 


r  dua 

E  h3 

r  4 

9  u 

2d2u 

+  s  , 

,r  + 

r 

4-ii 

be  +  ur. 

12b^ 

9S4 

n  1  U 

90 2  rj 

(7-4) 


The  quantity  a  is  the  surface  density  of  the  dome,  and  and 
?2  are  the  acoustic  pressures  in  regions  1  and  2,  respectively. 

The  surface  density  of  the  dome,  a,  must  be  expressed 
in  a  manner  that  accounts  for  surface  density  discontinuities  at 
structure  locations.  The  surface  density  of  the  dome  can  be 
written 

a  =  psh  +  ,  (7-5) 


where  og  is  density  of  dome,  h  is  thickness  of  dome,  M  is  mass 
per  unit  length  of  structure  element,  t  is  thickness  of  struc¬ 
ture,  A  is  cross-sectional  area  of  structure  element,  and 


6(9) 


0  >  (9  .  +  §) 
-  §  <  9  <  (9  .  +  §) 


9  <  (0  j  -  0 


(7-6) 
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One  additional  assumption  is  that  the  structure  element  does  not 
affect  the  tangential  D'Alembert  force. 

The  radial  and  tangential  velocity  components  of  the 
ftu  *u 

dome  are  given  by  W  =  =  -iru^,  and  V  =  ^-jr-  =  ixu^  ,  respec¬ 

tively,  for  harmonic  time  dependence.  When  the  above  are  sub¬ 
stituted  into  Eqs .  (7-3)  and  (7-4),  the  results  are 


x 


2 


hV  + 

L? 


Esh 


.2,, 

0  V 


(7-7) 


6  (9) 


tW  4 


Esh- 


5V 

5T 


"] 


5  vJ 


=  -i'uCPj  -  P2]  .  (7-8) 

The  velocity  components,  V  and  W,  are  the  total  velocity 
components  of  the  dome  as  are  the  pressure  fields  P^  and  P2  . 
Assume  that  the  total  velocity  components  and  pressure  fields 
may  be  written  in  two  parts: 


V  =  v  +  v  , 

(7-9) 

W  =  w  +  w  , 

(7-10) 

pi  =  p:  +  pi  , 

(7-11) 

p2  =  p2  +  p2  • 

(7-12) 

The  unbarred  quantities  are  the  non-structure  components,  and  the 
barred  quantities  are  the  components  due  to  the  presence  of  struc¬ 
ture.  If  Eqs.  (7-9)  -  (7-12)  are  substituted  into  Eqs.  (7-7)  and 
(7-8),  and  if  the  non-structure  solution  to  the  problem  is  sub¬ 
tracted  from  these  results,  the  resulting  equations  are, 
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2  ,  -  ,  Esh  0  Sw  n 

■"  ?ch  v  +  ~T  - 7  +  Ta  =  0  , 

s  b2  as2  d9 


(7-13) 


E  h  | 

-'J)2Dsh  W  -  'P2  ^  6  (c)t  (w  +  w)  +  — TJ-  +  W  1  + 


Esh  fa4w  ,  2S2w  ,  .1 


•  — Z  +  - 2  +  W 

l  Sc 

L  J 


+  w  I  =  -  ia>  T 


«■’  (Px  -  P23  •  (7-14) 


For  an  even  problem,  the  quantities  w,  v,  v  ,  and  p^. 
and  p.j  may  be  expressed  by  the  following  Fourier  series: 


v  =  ),D  srn  n9  , 
n  » 

n=o 


(7-15) 


CP 

=  S  Gn  cos  n0  , 


(7-16) 


Pi  =  2'-AnJn(kr)  +  EnNn(kr)3  cos  n6  ,  (7-17) 


P2  =  ^CnH]lp(kr)  cos  n9  . 


(7-18) 


When  Eqs .  (7-15)  -  (7-18)  are  substituted  into  Eqs.  (7-13)  and 
(7-14),  and  the  orthogonality  properties  of  the  sine  functions 
are  employed,  the  results  are 


2  ^  o  E  n  _ 

11  P  <,  ~  —j  n  D  -  — jy  G  =  0, 

b2  J  n  b2  n 


(7-19) 
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-I  .  h  G 
s  n 


2 - M-  2t  '  £ (w  +  w) 


cos  'fd-  +  — *■  : n  D  +  G  ]  — 
b2  n  n  "n 


F  h 

s‘  2-r_4 :  o_2; 


12b  n  n  n 

(7-20) 


u(D 


Equation  (7-16),  along  with  the  fact  that  the  structure  angular 
half-width  (?)  may  be  taken  as  small  as  desired,  allows  the  in¬ 
tegral,  I,  in  Eq.  (7-20)  to  be  evaluated  with  the  following  re¬ 
sults  : 


.  M,2  >aX  1  , 

1 - 5“  rr  (w  +  w)Q_a  cos  n6.  sin  nc 

j-1  J 


(7-21) 


Equation  (7-21)  when  substituted  in  Eq .  (7-20)  yields 

-  *  VVtb»  -  11  Wkb>  +  t*  CnH^>(kb)  -  -V  Dn  + 

b 


h» 


F  9  jmax  „ 

s  1  ,  h^  /4  o  2.-.VTIP  ,  Muu4  1  1  /  , ,  -  s  n 

71  1  +  7T7  'n  ~  2n  +1)]^Gn  +  7T(w+w)0=0  x 

b  12b  j=1  '  J 

(7-22) 


cos  n8 .  sin  n?  =  0  . 

J 

Boundary  conditions  on  the  surface  of  the  transducer 
must  be  matched.  It  is  important  to  realize  that  the  velocity 
prescription  of  the  transducer  face  has  not  been  changed  by  the 
addition  of  the  rib  structure  to  the  dome.  The  boundary  condi¬ 
tion  at  the  transducer  face  is 


(7-23) 
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(P]l  +  PX)  j 


(7-24) 


When  the  non-structure  result  is  subtracted  from  Eq .  (7-24),  the 
condition  to  be  met  is  that 


0  , 


(7-25) 


Anj;(ka)  +  Bnr(ka)  =  0 


(7-26) 


The  velocity  boundary  conditions  at  the  dome-fluid  interface  are 


=  ia'pW 

r=b  r=b 


(7-27) 


=  iuioW 

r=b  r=b 


(7-28) 


When  Eqs.  (7-15)  -  (7-18)  are  substituted  in  Eqs.  (7-27)  and 
(7-28),  the  results  of  matching  boundary  conditions  at  the  dome- 
fluid  interface  yield  the  following  relations, 


Vn(kb)  +  ®nNn(kb)  "  icP  5n  -  0  , 


=nHn:) ’ <kb>  -  5n  =  0 


(7-29) 


(7-30) 


In  Eq.  (7-16)  it  was  assumed  that  w  could  be  written 
as  a  Fourier  series.  Equations  (7-19),  (7-26),  (7-29),  and 
(7-30)  allow  Eq.  (7-22)  to  be  written  in  terms  of  only  one  unknown 
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coefficient,  namely  G^.  The  form  of  G  in  Eq  .  (7-16)  may  now  be 
specified  as 


G 


n 


n 


o  jmax 
Mx  2  V 

~TT”  L, 

j'l 


J_  (w  +  w)^=a  cos  n?. 

nl  “j  J 


sin  n? 


(7-31) 


where 


♦  = 
n 


-1 


>'cp!  Jn(kb)N;(ka)-Nn(kb)J^(ka)J  xcpHr(kb) - 

+  h 


u)2p 


-f  [1  +  -^-T(n4-2n2+l)l 

F~0 - ’ 

Es  n 

b1  12b2  '  ' 

b^(-x2psb2-Esn2) 

If  Eq.  (7-31)  is  substituted  in  Eq .  (7-16)  and  w  is  evaluated  at 
angles  a  =  c^,  then  w^  i<; 


x  jmax 

=  E'  £  W  Z  *n  sin  n?  E  <w  +  . 

K  n-o  '  J=1  J 


cos  ne .  cos  n8 


k  ’ 
(7-32) 


where  Z  is  the  "structure  coefficient"  defined  by  Z  =  > 


•  jk  2-t 


jmax 

V 

L-J 

("a 

j  =  l 

z 

2'rf 

CO 

E 

jk 


E  — — —  sin  nr  cos  n8  .  cos 
n  '  J 


n8 


k  ’ 


(7-33) 


(7-34) 


n=o 


98 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  78721 


jmax  jmax 

w,  =  V'  w  w.  v  . , 

;k  f-\  Jk  f-1.  M  Jk 

3= 1  J  J=1  J 


jmax 

w,  =  /  6  w„ 

%  ^  Jk  8, 

k  J  =  1 


6ik  = 


0  j  4  k 
1  j  =  k 


The  equation  for  w  may  be  expressed  in  matrix  form 


(1  "  ;11)  21 


Jmax  ,  1 


I  u 

2 

Jmax  * 


.  1  ,kmax  ~T2,kmax 


(1  -  U  ki-ax^ 

Vax-^l  Vax 


jmax 

*3  £  w3  , 

1  j  =  l  3j  J’1 

jmax 

:  w0  y  w0  C  .  2 

!  2  _  j  J  ’ 

.  .  .  ■  j=l„ 

„  -jmax 

z'  w0  t  •;  Ifmax 

kmax  A  'j  J)kmaX 


(7-35) 

and  the  matrix  may  be  inverted  by  standard  inversion  techniques. 

With  the  velocity  components  w  known,  the  coefficient  G  may 

9k  n 

be  determined  by  Eq .  (7-31)  and,  in  turn,  all  other  unknown  co¬ 
efficients,  D  ,  A  ,  B„ ,  and  C  can  be  determined.  This  allows 
computation  of  the  pressure  fields  ir  regions  1  and  2,  thus 
completing  the  solution. 
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8.  INFINITE  PLANAR  MODEL 

The  infinite  planar  dome- transducer  model  is  shown  in 
Fig.  8-1,  along  with  the  coordinate  geometry.  A  thin  plate  is 
located  a  distance  d  from  the  planar  "transducer."  The  trans¬ 
ducer  surface  has  a  continuous  normal  velocity  distribution  which 
simulates  the  particle  velocity  of  a  plane  wave  traveling  in 
a  direction  making  an  angle  ?  with  the  normal  to  the  plane. 

This  velocity  distribution  is  given  by 


v  .  v0elk  x  sin=  ,  (8-1) 

where  VQ  is  just  a  velocity  amplitude  factor.  The  pressure  field 
between  the  dome  and  transducer  will  be  denoted  by  p^,  and  that 
outside  the  dome  by  P2«  In  order  for  the  fluid  particle  velocity 
to  equal  the  transducer  velocity  at  y  =  0, 


_L_  |£|  =  v. 

HL-p  cyj 

y=0 


(8-2) 


The  plate  is  assumed  to  be  thin  enough  that  the  boundary 
conditions  at  its  surfaces  may  be  matched  at  the  middle  surface. 
One  condition  is  that  the  fluid  particle  velocity  at  the  inner 
and  outer  plate  surfaces  equal  the  plate  velocity.  If  the  plate 
displacement  from  equilibrium  is  denoted  by  u,  the  equality  of 
velocities  is  given  by 


1TC'  y=d 


1  rdp2n  OU 

1TP  y=d  Bt  ’ 


(8-3) 
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The  dynamic  behavior  of  the  plate  is  described  by  the  classical 
plate  equation.  In  this  case  the  driving  force  is  the  pressure 
difference  across  the  plate,  so  the  equation  takes  the  form 


PL  “  P2J 


y=d 


v  2 

h4  + 


Eh' 


.4 
o  u 


I2(l-v  )  3 x 


(8-4) 


Here  cc  is  the  plate  mass  density,  h  is  its  thickness,  E  is  Young's 
modulus,  and  v  is  Poisson's  ratio. 

The  solutions  to  the  wave  equation  for  the  pressure 
fields  inside  and  outside  of  the  dome  will  have  the  forms 


pl 


p2 


=  lA  eiTiy  +  B  e-^yieiKx  -2  •  -2  «-2 


r> 


K  =  k‘ 


=  C  e^y  eiKx  . 


(8-5) 

(8-6) 


Substituting  these  equations  into  the  boundary  condition  equations 
yields  expressions  for  the  unknown  coefficients  in  Eqs.  (8-4)  and 
(8-5).  The  final  results  for  the  pressure  fields  are 


Pl  *  'eiTl(7-d>  +  +  HK4) 

X  0 

Po  =  ^2  eiT1(y-d)  e1K>;  . 


cos  P (y-d) j elKx  (8-7) 


(8-8) 


where  r  =  c  h  ,  H  = 

C  * 


Eh' 


12(l-v  ) 


j  ,  K  =  k  sin  9  ,  ri  =  k  cos 


D  =  i^ie  *  ^  (-c  i  ^  4-  HK^)  sin  pd]  . 


.v  p 


It  is  seen  that  p^  is  the  sum  of  the  traveling  wave  which  is  trans¬ 
mitted  through  the  dome  and  a  standing  wave  which  always  has  an 
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antinode  at  the  plate.  The  average  intensity  passing  through 
the  plate  is  just 


I 


o 


(8-9) 


For  many  cases  of  interest  the  Hk  term  in  Eq.  (8-8)  is 
2 

small  relative  to  the  cn1  ,  and,  in  this  case 


D  ~  i- 


i-d 


sin  nd* . 

J 


Defining  R  as  the  ratio  of  D  in  Eq.  (8-9)  to  the  value  of  I  would 
have  with  no  dome,  h  =  0,  gives 

C  D 

R  =  1/  1  +  nh  —  sin  nd  [r.h  —  sin  nd  -  2  cos  nd*!,-.  (8-10) 

<-  p  D  J 

It  is  seen  that  R  is  a  periodic  function  of  d,  with  period 
~/k  cosS ,  and  that  its  amplitude  of  oscillation  decreases  as 
=  increases.  Starting  from  d  =  0,  R  increases  from  unity  until 
its  first  maximum  occurs  at  a  value  of  d  for  which 

tan  2  nd  =  2  rh  ps/c  .  (8-11) 

r\ 

Typically,  ps/o  =»  7.9  and  kh  «  1.8  x  10"  ,  so  for  9=0,  tan  2  -nd  = 
0.284.  For  this  small  value,  tan  2  nd  ~  0.14.  This  is  a  plate- 
transducer  spacing  of  0.023  wavelengths.  Maxima  occur  at  this 
value  of  d  and  at  every  half  wavelength  thereafter. 
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